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1.0  INTRODUCTION 

During  the  past  few  years  many  different  analysis  techniques  for 
turbulent  boundary-layer  flows  have  been  developed  and  validated.  While 
the  computational  aspects  of  these  developments  must  not  be  disregarded, 
the  real  key  to  the  understanding  of  how  to  properly  analyze  complex 
turbulent  boundary- layer  flows  Involves  advances  in  turbulent  transport 
modeling.  Reviews  of  many  of  the  various  proposed  mathematical  models 
of  turbulence  may  be  found  in  Refs.  1 through  6,  and  the  recent  books  by 
White  (Ref.  7)  and  Cebecl  and  Smith  (Ref.  8)  give  overall  coverage  of 
turbulent  boundary-layer  analysis  methods,  both  old  and  new. 

This  report  documents  another  turbulent  boundary-layer  analysis 
technique.  The  extended  mixing- length  hypothesis  is  first  developed; 
then  a baseline  turbulence  model  is  explored  and  the  baseline  model  is 
modified  to  account  for  effects  such  as  pressure  gradients  and  roughness. 
The  method  is  assessed  by  comparing  computed  results  with  experimental 
data  for  many  different  two-dimensional  and  axlsymmetrlc  turbulent 
boundary-layer  flows,  both  incompressible  and  compressible. 

2.0  ANALYTICAL  CONSIDERATIONS 

Many  recently  developed  computational  methods  for  turbulent  boundary- 
layer  flows  utilize  the  "eddy  viscosity"  concept  of  Boussinesq  and  apply 
Prandtl's  mixing-length  hypothesis  to  correlate  the  turbulent  shear 
stress  as  a function  of  the  local  mean  flow  field  in  the  boundary 
layer.  The  recent  compendium  by  Bushnell,  Cary,  and  Harris  (Ref.  9) 
provides  a list  of  such  programs.  The  method  of  Fatankar  and  Spalding 
(Ref.  10),  one  of  the  better  known  and  more  widely  documented  codes, 
correlates  the  mixing  length  in  terms  of  the  local  boundary- layer  thickness 
with  the  correlation  function  as  determined  from  flat-plate  flows  assumed 
applicable  to  general  flow  situations.  One  drawback  is  that  the  flat- 

I 

plate  correlation  represents  a reasonable  approximation  only  for  flows 
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with  moderate  pressure  gradients  and  Is  not  necessarily  applicable  for 
flows  with  large  pressure  gradients,  either  favorable  or  adverse. 

Drawing  a direct  relationship  between  the  mixing  length  and  the  local 
mean  flow  preempts  explicit  consideration  of  the  development  of  the 
turbulent  field,  and  such  methods  are  not  very  successful  In  predicting 
highly  "nonequlllbriiun"  flows  such  as  relamlnarlzatlon. 

The  present  analysis  utilizes  an  Integral  form  of  the  klnetlc- 
^nergy-of-^urbulence  (IKET)  equation  in  which  the  mixing  length  is  no 
longer  directly  correlated  to  the  local  flow  parameters  but  is  instead 
calculated.  The  additional  equation  thus  allows  the  "history"  of  the 
turbulent  state  to  be  considered  explicitly  and  the  mixing- length 
correlation  to  vary  as  the  turbulent  boundary  layer  develops.  This  so 
called  extended  mixing-length  approach  is  not  new,  having  been  developed 
by  McDonald  and  Camarata  (Ref.  11)  and  applied  by  McDonald  and  Fish 
(Ref.  12),  Shararoth  and  McDonald  (Ref.  13),  McDonald  and  Kreskovsky 
(Ref.  14),  and  Chan  (Refs.  15  through  17).  The  present  work  represents 
an  extension  of  the  same  approach  with  emphasis  on  compressible  flows. 
This  method  Is  similar  to  the  approach  of  Bradshaw,  Ferriss,  and  Atwell 
(Ref.  18)  with  the  exception  that  use  of  the  IKET  equation  eliminates 
the  pres sure- velocity  diffusion  process  and  necessitates  an  assumption 
as  to  the  distribution  of  mixing  length  across  the  boundary  layer. 

The  IKET  method  as  presented  in  this  report  lies  between  the  con- 
ventional mixing  length  and  the  one-equation  hydrodynamic  model  of 
turbulence,  to  use  the  hierarchy  of  Launder  and  Spalding  (Ref.  1). 

2.1  GOVERNING  EQUATIONS 

The  present  analysis  employs  the  compressible,  time-averaged 
boundary- layer  equations  for  two-dimensional  or  axlsymmetrlc  geometries 
as  derived  by  Vaglio-Laurln  (Ref.  19).  As  shown  in  the  sketch,  the  body 
surface  is  defined  by  y = 0,  and  x is  defined  as  along  the  body  surface. 
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For  axlsymmeCrlc  bodies,  r(x)  is  the  radius  at  x and  y ° 0.  The  velocity 
components  are  taken  to  be  u and  v in  the  x-  and  y-  directions,  respec- 
tively. 


The  governing  equations  of  motion  are  as  follows  (see  Nomenclature  for 
terminology) : 


*An  additional  term  for  wall  roughness,  explained  In  Section  2.4. 
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The  usual  expressions  for  the  mean  and  fluctuating  parts  of  the  depen- 
dent variables  are  used;  e,g., 

u = u + u*"  (6) 

Implicit  In  Vagllo-Laurln's  derivation  of  the  above  equations  are  the 
following  stipulations: 

a.  The  rates  of  change  of  the  mean  flow  properties  In  the  x 
direction  [0(1)]  are  smaller  than  the  rates  of  change  In 
the  y direction  [0(6  ^)]  by  an  order  of  magnitude. 


b.  Mean  squares  and  products  of  the  turbulent  fluctuations  are 
0(6);  that  Is,  the  turbulent  level  Is  small.  The  terms 
Involving  mean  squares  of  the  velocity  fluctuations  are 
taken  to  be  negligible,  which  Is  valid  for  high  Reynolds 
number  flows  with  zero  or  favorable  pressure  gradients. 


c.  The  time-averaged  molecular  transport  quantities  are  approxi- 
mated by  those  pertaining  to  the  mean  flow  properties;  todeed, 
even  the  latter  are  negligible,  except  very  near  the  wall, 
when  compared  with  terns  Involving  the  turbulent  transport 
quantities. 


If  subscript  w denotes  the  wall  and  subscript  e denotes  the  cuter 
edge  of  the  boundary  layer,  the  associated  boundary  conditions  on  the 
aforementioned  equations  are 

- .V  - . C ' ‘ 

y = 0:  u = u'v'  = p'v'  = v'h'  = 0 


h = h^  (prescribed  wall  enthalphy) 

~ = 0 (adiabatic  wall) 
iy 

V = 0;  V p =■  (v  p)  (transpiration) 
w 
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llm  y -»•  »:  = u^v'*  = p'v'  = 0 


u 

h 


e 

which  reflects  the  no-sllp  condition  at  the  wall  as  well  as  a prescribed 
wall  enChalphy  (or  temperature)  or  an  adiabatic  wall.  The  normal 
momentum  equation  Eq.  (3}j  implies  that  the  static  pressure  variation 
across  the  boundary  layer  is  negligible  and  that  the  static  pressure 
distribution,  p(x),  is  an  external  imput  to  the  bound ary- layer  analysis 
from  a separate  Inviscld  analysis  or  from  experimental  data.  The  outer- 
edge  velocity,  U^,  as  well  as  the  outer-edge  enthalpy,  h^,  must  be 
detemlned  from  an  inviscld  analysis  consistent  with  the  Imposed  pres- 
sure distribution.  The  gas  model  adopted  for  the  present  study  is  ther- 
mally and  calorlcally  perfect  air  (or  nitrogen)  obeying  the  relations 


p ■=  p R T 


(9) 


h = c T 
P 


(10) 


with  the  specific  heat  ratio  y = 1.40,  the  gas  constant  R = 1,716  ft- 
lbf/siug-‘’R  for  air  (1,776.29  for  nitrogen),  and  the  specific  heat  at 
constant  pressure  = 6,006  f t-lbf /slug-®R  for  air  (6,217  for  nitrogen). 
The  laminar  viscosity  is  taken  to  obey  Sutherland's  law, 

,-8  T Ibf-sec 


2.27  X 10 


T + 198.6 


ft^ 


(11) 


(T  is  in  degrees  Rankine).  The  laminar  Prandtl  number,  Fr,  is  assumed  to 
have  a constant  value < of  0.71  across  the  boundary  layer;  l.e. , 

Pr  “ 0.71>  (12) 


Following  the  development  given  in  Appendix  A of  the  report  by 
Laster  (Ref.  20),  the  continuity,,  energy,  and  streamwlse  momentum 
equations  can  be  combined  to  yield  the  compressible  time-averaged 
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turbulent  kinetic  energy  equation  for  two-dimensional  .boundary- layer 
flows.  The  approach  of  Laster  In  conjuctlon  with  equations  equivalent 
to  Eqs.  (1),  (2),  and  (4)  yields  the  compressible,  time-averaged  tur- 
bulent kinetic  energy  equation  valid  for  two-dimensional  or  axlsymmetric 
flows,  which  may  be  written  in  the  form 

1/r^  3/3x  (\'I2  r^  ^ vrq^^+  3/3y  |l/2  (Fv  + p"  v")  q^ 

= - (p  u''  v'  + p"  u''  v"*)  3u/3y 
+ p**  3v‘*/3y  - p"*  u"  (u  3u/3x  + v 3u/3y)  - p (u**)^.  3u/3x  - p (v")^  3v/3y 


+ 3/3y  (\/2  p'  v^  q^jl  - 3/3y  ^p^v' 


+ 1/2  p v"  q^  +1/2  p 


' v') 


k 8/Sy  [l 


+ 3/3y  3/3y  |l/2q2  + (vO^ 


vO^]  I - P t 


(13) 


where 


P E 


' m ' 


(14) 


represents  the  viscous  dissipation  of  turbulent  energy  in  conventional 

2 ' 

tensor  notation.  Above,  1/2  q denotes  the  turbulent  kinetic  energy,  k, 
defined  as  ’ i : 

k = 1/2  q^  ■=>  1/2  [(u^)^  + (v')^  + (wO^] 


A discussion  of  the  physical  meaning  of  each  term  in  Eq.'  (13)  Is  given  ' 
in  section  IV  of  chapter  II  In  the  report  by  Laster  (Ref.  20)  which  also 
contains  a derivation  of  the  compressible  turbulent  kinetic  energy 
equation  In  tensor  notation.  The  paper  by  Bradshaw  and  Ferrlss  (Ref. 

21)  and  the  book  by  Cebecl  and  Smith  (Ref.  8)  contain  similar  deriva- 
tions and  discussions.  Equation  (13)  has  been  derived  using  no  assump- 
tions or  restrictions  other. than  the  classical  boundary— layer  order— of— 
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magnitude  analysis  used  to  eliminate  certain  terms  from  the  complete 
turbulent  kinetic  energy  equation.  The  turbulent  kinetic  energy  equa- 
tion [Eq.  (13)]  can  be  formally  Integrated  from  the  wall  (y  ^ 0)  to  the 
edge  of  the  boundary  layer  (y  = 6)  to  obtain 

o ^ ^/2  r^  p u q^^  dy  + [i/2  (p  'v'+  p‘*v'‘)  <1^ 


>6 


- /°  (p  u^v"  + p'  u'  v')  3u/9y  dy  + /'"  p'  3v'/ 9y  dy 
o o 


- r 

0 


^p'u'  (u  9u/9x  + V 9u/9y)  + p(u‘')^  9u/9x 


+ P (y')^  9v/9yJ  dy  + (l/2  p**  v"*  q^)^ 


- (p"v^+  1/2  p 


v'  q'^  + 1/2  p 


q^v') 


+ {y  3/3y  [t/2  q^  + (v-)^]}o  - / 


6 _ 
p e dy 


(16) 


Leibnitz's  rule  (Ref.  22)  for  differentiation  of  an  Integral  depending 
upon  a parameter  must  be  used  since  the  boundary-layer  thickness,  5,  Is 
a function  of  the  streamwlse  distance  x;  l.e.,  5 = 5(x).  Thus  the  first 
Integral  In  Eq.  (16)  can  be  written 

9/9x  ^1/2  r^  p u q^^  dy  = d/dx  r^  1/2;  p u q^  dy 


- 1/2  r^  Pg'Ugq^  d6/dx 


(17) 


where  the  subscript  e denotes  conditions  at  the  edge  of  the  boundary 

layer  where  y =6,  Since  all  turbulent  quantities  such  as  q^,  p'v't 
etc.,  are  Identically  zero  at  the  wall  (due  to  viscous  damping)  and  ' ' 

since  9/9y  q^  + (v'')^Jls  Identically  zero  at  1»th  the  wall  and  the 
boundary-layer  edge,  Eq.  (16)  may  be  written  as 
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l/r^  d/dx  / 1/2  p u q^dy  = - / (p  u'v'  + p'u'v') 

0 0 . ' 

(S  <S  / f 

3u/3y  dy  + / p'av'/Sy  dy  - / ) p'u'fu  9u/3x  + 

0 o' 

V 3ll/3y]  + ^(u')^  9u/3x  + p (v')^  3v/3y  dy  - 
6 _ 

/ p e dy  + E (13) 

with  the  term  E defined  as 

E = ^1/2  q2  ("^  u dd/dx  - p v)  - p'v'-  1/2  "p  v''q^  , 

- 1/2  p'  q2v'  1 (19) 

The  subscript  e denotes  that  all  quantities  in  the  term  E are  to  be 
evaluated  at  the  edge  of  the  boundary  layer.  Equation  (18)  Is  the  IKET 
(Integral  kinetic-energy-of-turbulence)  equation.  A discussion  of  the 
physical  meaning  of  the  various  terms  in  Eq,  (19)  will  be  given  in 
Section  2.5. 

2.2  TURBULENCE  MODEL 

The  governing  equations  [Eqs,  (1),  (2),  (4),  and  (18)]  presented 
and  examined  in  the  preceding  section  do  not  in  their  present  form  con- 
stitute a complete  system,  since  there  are  more  dependent  variables  than 
equations.  In  particular,  the  fluctuating  quantities  Introduced  by  the 
time-averaging  process  must  be  functionally  related  to  the  mean  flow 
variables  or  additional  partial  differential  equations  provided  for  the 
fluctuating  quantities;  i.e.,  the  problem  of  closure  must  be  faced. 

Cebeci  and  Bradshaw  (Ref.  23)  discuss  In  an  admirable  manner  the  problem 
of  closure  and  Its  ramifications.  Turbulence  models,  plausible  simpli- 
fying assumptions  involving  mathematical  modeling  of  various  turbulence 
quantities,  are  used  to  achieve  closure  and  hence  to  reduce  the  number 
of  dependent  variables  to  the  number  of  equations.  It  is  via  this 
modeling  process  (content  and  completeness)  that  this  analysis  differs 
markedly  from  the  previous  extended  mixing- length  Investigations  reported 
(Refs.  11-17). 
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The  turbulent  shear  stress,  for  compressible  flow  is  given  by 

Tj.  = - p - p'u'v"  (20) 

Following  the  assumptions  and  approximations  of  Bradshaw  and  Ferrlss 
(Ref.  21),  namely  that  p"‘u'’v‘'  <<  p u^v",  It  Is  assumed  that 

- p u^v"  (21) 

The  Kolmogorov-Frandtl  model  of  turbulence  as  presented  by  Wolfshteln 
(Ref.  24)  Is  adopted  so  that 

= Pj  3u/3y  (22) 


with  the  turbulent  viscosity  defined  by 

p^  = Cy  P 


(23) 


where  is  an  empirical  constant,  k is  the  kinetic  energy  [defined  by 
Eq.  (15)3  ' length  scale  for  turbulent  shear  stress. 

The  dissipation,  D,  of  turbulent  energy  Is,  following  Wolfshteln  (Ref. 
24), 


D = p e 


Cjj  p k 


3/2 


(24) 


where  Is  an  empirical  constant  and  Is  the  length  scale  for  dissi- 
pation. In  accord  with  Bradshaw,  Ferrlss,  and  Atwell  (Ref.  18)  and 
Townsend  (Ref.  25), 


= aj  - 2ai  k = x^/p 

(25) 

(u')^  = 32 

(26) 

(v‘')2  = a3  q2 

(27) 

~~2 

« (1  - a,  - a.)  q 

(28) 
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The  numerical  values  of  a^,  a2>  and  a^  are  generally  ascribed  constant 
values  which  are  Independent  of  the  x-  and  y-coordlnate  directions.  The 
work  of  Rose  and  Murphy  (Ref.  26),  however,  established  that  a^  is  not  a 
constant  but  Is  a function  of  the  distance  normal  to  the  wall  (the  y- 
coordlnate  in  this  report);  hence,  a^  = a-|(y)*  , \ 

Terms  involving  the  fluctuating  density,  p',  are  modeled  following 
Bradshaw  and  Ferris s (Ref.  21)  as 

p '/p  R;  uu^/h 

so  that  ' 

r a,  2a-  "p  uk 

p'u'  (pu/h)  (u"")^  = = 

li  h 

and  ^ 

' p U U^v'  - T U 

p'v'  % — IT = — — 

h h 

The  fundamental  assumptions  leading  to  Eq.  (29)  are  that  (1)  the  fluctua- 
tion of  total  enthalphy  is  much  less  than  the  fluctuation  of  static 
enthalpy  and  (2)  the  pressure  fluctuations  are  small  If  the  Mach  number 

fluctuation  Is  much  less  than  unity.  It  Is  also  assumed  that  the  term 
2 

(y  - 1)M  is  no  greater  than  unity  where  M is  the  local  .time-averaged 
Mach  number.  Shamroth  and  McDonald  (Ref.  13)  discuss  the  relative 
Importance  of  the  pressure  dilatation  term,  p"  Sv^/dy,.  and  conclude  that 
it  is  small  relative  to  the  other  terms  In  the  turbulent  kinetic  energy 
equation.  Combination  of  all  the  above  assumptions  essentially  re- 
stricts the  formal  applicability  of  the  resulting  analysis  to  supersonic 
and  low  hypersonic  flows  with  free-stream  Mach  numbers  less  than  ap-  v 
proximately  five. 


(29) 

1 

(30) 


: (31) 
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The  turbulent  heat  flux,  p v'^h",  Is  modeled  via  use  of  a turbulent 

Prandtl  Number,  Pr^  (baaed  on  the  use  of  static  enthalphy),  by  which  the 

turbulent  thermal  conductivity  Is  related  to  the  turbulent  viscosity  as 

Pr  = c The  turbulent  Prandtl  number,  Pr.  , is  taken  to  be  a 

t p t t t 

constant  value  of  0.90  across  the  entire  boundary  layer. 


The  time-averaged  quantities  (such  as  u'v^)  appearing  In  the  IKET 
equation  [Eq.  (18)3  related  to  the  mean  flow  quantities  and 

the  various  length  and  structural  scales  and  constants  introduced  by  the 
turbulence  model.  In  particular,  Eqs,  (21),  (22),  (23),  and  (25)  can  be 


combined  to  yield 

- u"*  v"*  = Tj,/p  = ^y^/p^3u/ay  = 3u/3y  = 2a^k  (32) 

which  results  In 

k^^^  = (c^/2aj  3u/3y  (33) 

so  that  the  turbulent  viscosity,  y^,  given  by  Eq.  (23)  becomes 

y^  =(A2a^)F  (34) 

and  hence  the  turbulent  shear  stress  may  be  written  as 

Tt  “ yt  3u/3y  = ^cj/2a^^7  (3u/3y)^  (35) 


The  dissipation  of  turbulent  energy,  given  by  Eq.  (24),  can  now  be 
expressed  as 

D « e = (V^V  ^ 
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Using  the  aforementioned  relations,  one  can  combine  terms  in  Eq, 
(18),  the  IKET  equation,  to  produce 

1/rJ  d/dx  f rJ  ? u {9u/3y)^  dy 


- / ^ [‘  - <V«D>] 

- / p(£y/2ai)^  Ou/3y)^  [232  (1  + u^/h)  au/9x 


+ 2a2(u  v/h)  3u/3y  + 233  3v/3yJ  dy  + E/C^  + E^ 


(37) 


The  additional  term  involving  E^  (reflecting  wall  roughness  effects) 
will  be  explained  in  Section  2.4.  A discussion  of  the  physical  meanings 
of  the  terms  E and  E^  will  be  given  in  Section  2.5. 

2.3  LENGTH  SCALES,  STRUCTURAL  SCALES.  AND 
EMPIRICAL  CONSTANTS 


The  implementation  and  effective  use  of  Eq,  (37)  is  dependent  upon 
the  length  scales  (£^  and  £jj) , the  structural  scales  (a^ , and  a^), 
and  empirical  constants  (Cj^  and  C^).  A two-layer  model  of  the  turbulent 
boundary  layer  is  adopted  following  the  classical  inner-outer  region 
approach  in  which  separate  functional  relationships  are  prescribed  in 
each  region  with  continuity  of  the  functions  between  each  region. 


For  the  inner  region,  which  includes  the  viscous  sublayer,  the 
length  scales  and  the  structural  scale  are  taken  as 


= y Dp 

(38) 

= y Dd 

(39) 

^1  “ at  Da 

(40) 
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where  D , D_,  and  D are  damping  functions  which  account  for  viscous 
- p’  D a r o 

sublayer  effects  and  which  are  assumed,  following  Wolf shteln  (Ref.  24), 

to  be  of  the  form 

= 1 - exp  (-A^  R^) 

(41) 

Dj)  = 1 - exp  (-Ad  R(.) 

(42) 

= 1 - exp  (-Aa  R^) 

(43) 

with  empirical  constants  A , A-,  and  A to  be  determined. 

In  Eq.  (40) 

the  quantity  represents  the  fully  turbulent  value  of  the  structural 
scale  a^,  whereas  denotes  a local  turbulence  Reynolds  number  defined 
as 

-.1/2 

p k y 

R 

V 


The  length  scales  and  the  structural  scale  a^  for  the  outer  region 
are  of  the  form 


1 = Y 


ai  = at 


(45) 

(46) 

(47) 


where  Is  an  empirical  constant,  5 Is  the  boundary- layer  thickness, ' 
and  Y Is  the  outer  region  value  of  the  length  scale  Continuity  of 

the  functional  relationships  is  used  to  couple  the  inner-  and  outer- 
region  values  of  and  a^ . That  Is,  continuity  Is  maintained 

between  Eqs.  (38)  and  (45),  Eqs.  (39)  and  (46),  and  Eqs.  (40)  and  (47). 
The  addition  of  the  IKET  equation  allows  one  parameter  normally  taken  as 
constant  to, be  solved  for  In  the  streamwlse  (x)  direction.  Conventional 
mixing-length  hypothesis  turbulence  models  take  1^/6  to  be  a constant 
(usually  0.09  and  written  as  X/6).  The  extended  mixing-length  hypothe- 
sis turbulence  model  examined  In  this  report  allows  Y to  be  a variable 
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which  is  determined  (i.e.,  calculated)  by  the  IKET  equation.  The  outer- 
region  value  of  the  length  scale  £ , namely  Y,  is  the  one  parameter  to 
be  determined  from  Eq.  (37)  and,  hence,  controls  the  streamwise  develop- 
ment of  the  turbulent  shear  stress  in  such  a way  that  the  history  of  the 
turbulent  state  is  considered  explicitly. 

The  values  of  the  nine  empirical  constants  (C  , C„,  a , a_,  a-,  A , 

p D t 2 3 p 

Ajj,  A^,  and  A^)  introduced  previously  can  be  plausibly  deduced  using 
classical  turbulent  boundary-layer  theory  as  well  as  recent  experimental 
results  involving  compressible  turbulent  boundary-layer  flows.  Appendix 
A contains  the  details  of  such  a process.  The  deduced  values  of  the 
constants  together  with  the  values  used  by  Wolfshteln  (Ref.  24)  are 
given  in  Table  1. 


Table  1.  Turbulence  Model  Empirical  Constants 


Present  Work 

Wolfshteln® 

s 

0.2383 

0.220 

Cd 

0.3777 

0.416 

^t 

0.150 

a2 

0.566 

®3 

0.150 

0.016 

0.016 

0.18885 

0.263 

^a 

0.0469 

0.2069 

^ef.  24 

2.4  TURBULENCE  MODEL  MODIFICATIONS 

The  values  of  nine  empirical  constants  contained  in  Table  1 repre- 
sent a "baseline"  turbulence  model  and  should  be  viewed  as  valid  for 
zero  pressure  gradient,  no  transpiration,  and  smooth-wall  boundary-layer 
flows.  Because  of  the  intimate  relationship  which  exists  between  the 
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surface  condition  (roughness, . transpiration,  etc),  the  pressure  gradient, 
and  the  structure  of  the  turbulent  boundary  layer  near. the  surface,  some 
of  the  constants  listed  In  Table  1 can  be  expected  to  have  magnitudes  : 
different  from  the  baseline  case. 

As  utilized  in  Appendix  A the  hypothesis  of  Morkovin  (Ref,  27) 

(which  states  that  the  turbulent  motion  and  structure  in  a compressible 
boundary  layer  should  be  equivalent  to  those  of  an  Incompressible 
boundary  layer  as  long  as  the  turbulent  motion  Itself  Is  Incompressi- 
ble — usually  implying  a local  Mach  number  less  than  five)  permits  the 
employment  of  the  vast  incompressible  data  on  turbulent  boundary- layer 
structure  for  turbulent  compressible  boundary  layers.  Thus  the -effects 
of  roughness,  transpiration,  and  pressure  gradient  which  have  been 
extensively  and  carefully  ascertained  from  incompressible  experiments  •' 
can  be  used  to  Infer  changes  In  the  nine  constants  of  Table  1 from  their 
baseline  values. 

These  data  suggest,  at  least  within  the  context  of  the  mixing- 
length  concept,  that  the  structure  near  the  wall  Is  most  significantly 
affected.  This  corresponds  in  the  inner-outer  region  analysis  employed 
in  this  paper  to  the  Inner-region  structure  which  is  controlled  by  the 

, a^,  and  a^  are 

assumed  to  be  essentially  "universal"  constants  which  are  Invariant 
, over  a wide  range  of  conditions.  The  results  of  Rose  and  Murphy  (Ref. 

26)  for  a^  tend  to  confirm  this,  and  by  inference  the  Morkovin  (Ref.  27) 
hypothesis  supports  this  for  82  and  a^.  The  constants  and  apply 
across  the  entire  boundary  layer,  both  the  inner  and  outer  regions,  and 
hence  should  vary  only  slightly  and  in  particular  should  be  Independent 
of  the  values  of  the  inner-region  constants  A^,  A , and  A . The  remain- 
Ing  constant,  is  germane  only  in  the  outer. region  and  should  not  be 

directly  coupled  to  the  inner-region  constants.  . . 


constants 


and  A 


The  structural  constants  a^ 
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As  demonstrated  In  Appendix  A,  the  following  equations  may  be 
deduced  using  the  previously  defined  expressions  for  the  present  turbu- 
lence model  and  the  existing  data  base  for  turbulent  boundary- layer 
flows; 


^ = X/k 


(48) 

(49) 

(50) 


1/2 


(2a.) 2 


(51) 


^eff 
2a^  ic 


(A^  C ) 


3/2 


(52) 


0.016 


(53) 


Here  k Is  the  von  KArmAn  constant  and  A ..  Is  the  effective  value  of  the 

eft 

constant  In  the  classical  van  Driest  (Ref,  28)  damping  expression.  The 
term  A^^^  may  be  Interpreted  as  the  effective  thickness  of  the  viscous 
sublayer  expressed  In  terms  of  the  Inner-layer  coordinates.  It  Is  well 

. '..L:  i ' ‘ r 

known  (e.g.,  van  Driest)  that  effects  such  as  transpiration,  pressure 
gradient,  and  roughness  can  be  simulated  by  changing  the  damping  expres- 
sion since  these  effects  tend  to  alter  the  sublayer  thickness.  Kays' aiid 
Moffat  (Ref.  29)  and  Healzer,  Moffat,  and  Kays  (Ref.  30)  present  empiri- 
cal correlations  of  the  above  effects  expressed  In  terms  of  the 

classical  van  Driest  damping  parameter.  An  examination  of  Eqs.  (51) 
through  (53)  indicates  that  A^  reflects  a dependence  upon  A^^^  but  that 
Ajj  and  A^  do  not.  Since  the  structure  of  the  Inner  region  must  reflect 
the  effects  of  transpiration,  pressure  gradient,  and  roughness,  the 
expressions  A^^  and  A^  as  given  by  Eqs.  (51)  and  (53)  must  be  altered. 

The  forms  chosen  for  A^^  and  A^  will  be  given  later  in  this  report.  At 
this  point  the  realization  that  and  A^  as  well  as  A^  must  be  func- 
tions of  A^--  is  sufficient, 
eff 
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In  an  inner—  and  outer— region  model  of  the  turbulent  boundary  layer 
the  single  most  Important  parameter  is  the  effective  thickness  of  the 
viscous  sublayer,  Although  it  constitutes  only  a small  fraction 

of  the  boundary-layer  thickness,  the  sublayer  is  the  region  wherein  the 
major  change  in  velocity  occurs  and  where  (except  for  low  Prandtl  number 
fluids)  the  most  resistance  to  heat  transfer  takes  place.  Viewed  in  the 
context  of  the  viscous  sublayer,  experimental  evidence  supports  the 
concept  that  a favorable  gradient  (dp/dx  negative)  results  in  Increased 
sublayer  thickness,  an  adverse  pressure  gradient  (dip/dx  positive) 
results  in  decreased  sublayer  thickness,  transpiration  (blowing  or 
suction)  affects  (decreases  or  Increases)  the  sublayer  thickness,  and 
surface  roughness  results  in  decreased  sublayer  thickness. 

Kays  and  Moffat  (Ref.  29)  recommend  the  following  equations  for  the 
effective  equilibrium  value  of  the  viscous  sublayer  thickness  as  a 
function  of  pressure  gradient  and  wall  transpiration  velocity: 


A+  = _ 

26.0* 

%ff 

/ D+  \1 

a 

1 + CV+ 

+ 1 

L \ w/J 

where  v is  the  wall  transpiration  velocity  normalized  by  the  friction 

velocity  u^  =»  and  p is  the  dimensionless  pressure  gradient 

defined  as  constants  a and  b are  functions  of  p"*" 

and  V and  are 
w 

a = 7.1  for  v^  0 ' otherwise  a = 9.0 
b = 4.25  for  p^  £ 0 otherwise  b = 2.0 

c = 10.0  for  p^  ^ 0 otherwise  c =»  0.0 


*24.0  as  given  by  Ref,  29  but  taken  as  26.0  in  this  report  for 
consistency  with  other  empirical  values. 
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The  effects  of  pressure  gradient  and  transpiration  can  be  seen  In  Fig. 

1,  where  Eq,  (54)  is  parametrically  plotted.  Figure  1 graphically 

illustrates  the  high  values  which  A assumes  in  the  presence  of  strong 

ett  ^ 

favorable  pressure  gradients. , These  large  values  of  result  in  the 

viscous  sublayer  overwhelming  the  entire  boundary  layer  (relamlnariza- 
tlon)  and  in  decreased  heat  transfer  because  of  the  Increased  resistence 
within  the  sublayer.  Thus  a considerable  amount  of  physics  is  expressed 
by  Eq.  (54) . 


Figure  1.  A^ff  as  a function  of  pressure  gradient 
and  transpiration  [Eq.  (54)] . 


Equation  (54)  is  essentially  based  upon  the  assumption  of  inner-  .. 
region  equilibrium  and  as  such  does  not  describe  the  manner  in  which  ;; 

varies  while  changing  from  one  equilibrium  condition  to  another  or 
while  seeking  an  equilibrium  condition.  The  recent  paper  of  Horstman 
(Ref.  31)  evaluates  several  approaches  and  concludes  that  a lag  model  in 
which  A^j^  is  a function  of  weighed  upstream  values  of  p^  gives  results 
superior  to  those  of  a rate  model  in  which  the  rate  of  change  of  A^^^ 
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is  proportional  to  the  difference  between  and  the  equilibrium 

value.  In  the  lag  formulation,  p"^ (x)  in  Eq.  (54)  Is  replaced  by  p^(x) 
where 

P'^Cx)  “ / u(?)  d?  (55) 

x-2X£6 

and 

0,(5)  = —1—  exp  [-  U - (x  - X»6)|2/2a^]  (56) 

/lir  a 

with 

a - (57) 

where  X^  is  the  lag  length  parameter.  The  basic  effect  of  the  lag 
formulation  is  not  to  diminish  or  enhance  the  effect  of  pressure  gra- 
dient but  to  shift  the  effect  downstream.  Horstman  makes  no  attempt  to 
prescribe  X^  as  a function  of  x,  but  he  does  suggest  that  X^  should  be 
much  smaller  for  favorable  pressure  gradients  than  for  adverse  pressure 
gradients.  Values  of  20  and  10  were  used  for  adverse  and  favorable 
pressure  gradients,  respectively.  These  values  of  X^^^  do  not  represent 
the  optimum  physics  for  every  flow  condition,  but  they  do  represent 
reasonable  values.  Horstman  (Ref.  31)  states  that  additional  experi- 
mental data  are  needed  for  values  of  p*^  between  0,01  and  0.06  for  a wide 
range  of  Reynolds  numbers  in  order  to  develop  a realistic  correlation. 

The  results  of  Glowackl  and  Chi  (Ref,  32),  as  presented  by  Horstman, 
indicate  that  the  von  Karman  constant,  k.  Is  not  a universal  constant 
but  is  mildly  dependent  upon  the  pressure  gradient.  In  particular  they 
suggest  that 

K = 0.435  + 0.182257  [l  - exp (-0,320686)]  (58) 
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where  B is  a dimensionless  pressure  gradient  defined  as  5.  (dp/dx)/T 
with  5^  as  the  kinematic  displacement  thickness.  As  with  in  Eq. 

(54)  this  represents  an  equilibrium  value  of  ic.  The  lag  analysis  can 
be  applied  to  ic  by  replacing  @ In  Eq.  (58)  with  8 where 

I(x) 

with  aj(c)  defined  In  Eq.  (56).  Glowackl  and  Chi  recommend  that  Eq.  (58) 
be  used  for  flows  with  the  dimensionless  pressure  gradient  8 > 0 and 
imply  that  the  equation  can  be  used  for  6 < 0,  In  this  report  for  8 5.  0 
(dp/dx  <_  0.0)  the  von  Karman  constant  k was  taken  to  be  defined  by  Eq. 
(59)  with  the  stipulation  that  ic  must  he  greater  than  0.05. 

The  value  of  A^^^  (and  perhaps  ic)  calculated  in  the  above  manner 
must  now  be  related  to  the  Inner-region  empirical  constants  as  given  In 
Eqs.  (48)  through  (53). 


/ 


ft)(c)  8(C)  dc 


(59) 


x-2X^6 


Equations  (48)  through  (50)  and  Eq.  (52) , ref lect,  in  a logical 
manner,  the  importance  of  k and  0*1  the  constants  ,C^,  Cj^, 

and  A and  are  used  In  this  analysis.  However,  as  previously  mentioned, 

. d 

Eqs.  (51)  and  (53)  for  A^  and  A^  do  not  adequately  reflect  changes  in 

the  inner  region  because  of  transpiration,  pressure  gradient,  and 

roughness.  Hence,  changes  (as  predicted  by  the  lag  approach,  ' 

for  example)  alter  the  length  scale  for  mixing  and  must  be  reflected  in 

the  formulation  used.  Because  of  differences  in  the  behavior  of  Eq. 

(38)  and  the  classical  van  Driest  relationship  as  the  wall  is  approached, 

no  unambiguous  relationship  is  evident,  but  Inspection  of  the  two  models 

suggests  that  A and  should  be  related  in  an  inverse  manner.  The 

y err 

relationship  used  in  this  analysis  Is 


10.816 


(60) 
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which  reduces  to  A = 0.016  for  = 26.0.  Equation  (60)  has  been  ade- 

VI  eff  ^ 

quate  for  a wide  range  of  values  of  In  a similar  vein,  A^j,  the 

damping  constant  for  the  dissipation  length  scale,  has  been  taken  as 


These  equations  were  verified  by  computer  experimentation  and  were  used 
unaltered  for  a wide  variety  of  conditions  examined  In  Section  3.0  of 
this  report. 


The  modifications  made  to  the  baseline  turbulence  model  to  account 
for  the  effects  of  surface  roughness  are  more  extensive  than,  although 
similar  to,  those  required  for  pressure  gradient  effects.  The  roughness 
model  formulated  uses  the  approach  of  Healzer,  Moffat,  and  Kays  (Ref. 

30)  in  conjunction  with  the  analysis  of  Finson  (Ref.  33).  The  usual 
method  of  examining  and  presenting  boundary- layer  data  over  rough  sur- 
faces is  to  classify  the  surface  as  smooth,  transitional,  or  fully 
rough  according  to  the  local  value  of  the  roughness  Reynolds  number 


(62) 


where  u is  the  wall  friction  velocity.  The  surface  is  considered 
w 

smooth  aerodynamically  If  R is  less  than  5.0  and  the  surface  is  fully 
rough  for  R^  greater  than  55.0  (although  Schllchtlng  (Ref.  34),  for 
example,  suggests  70.0);  otherwise,  the  surface  is  transitionally  rough. 
For  boundary- layer  flows  with  R^  less  than  5.0  the  surface  is  taken  as 
smooth  with  no  resulting  effects  of  roughness. 


The  model  for  the  transitional  roughness  regime  will  now  be  devel- 
oped. Following  Healzer  et  at.  (Ref.  30),  the  effect  of  roughness  in 
this  regime  can  be  viewed  as  reducing  the  sublayer  thickness  which.  In 
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turn,  can  effectively  be  controlled  by  decreasing  hence 

[through  Eq.  (60)]  Increasing  A^.  The  functional  relationship  Is  taken 
as 

4.007  In  K (63) 

eff  * f 

where  ; ' 

f = 1.299  ("uniform"  roughness) 

2.061  (sand  grain  roughness) 

3.196  (commercial  roughness) 

+ + 

and  A In  Eq.  (63)  is  taken  as  from  Eq,  (54).  Equation  (60)  Is 

then  used  to  relate  A^^^  to  A^.  The  quantity  f in  Eq.  (63)  is  a func- 
tion of  the  uniformity  of  the  roughness.  The  value  of  f varies  from 
3.196  for  commercial  roughness,  which  Is  completely  random  In  shape  and 
size  distribution,  to  1.299  for  uniform  shape,  size,  and  packing  distri- 
bution. Physically,  this  corresponds  to  a more  rapid  thinning  of  the 
sublayer  by  the  more  random  roughness.  Figure  2 Illustrates  Eq.  (63) 
for  each  of  the  three  values  of  f as  a function  of  R with  A***  taken  as 

T 

26.0.  The  quantity  (4.007  - JlnR^)/f  cannot  exceed  1.00  in  value  as  this 
would  correspond  to  a thickening  of  the  viscous  sublayer  by  roughness,' 
which  Is  contrary  to  physical  reality. 

The  regime  of  fully  rough  turbulent  flow  (R^  > 55.0)  is  predicted 
using  the  approach  of  Finson  (Ref,  33).  Physically,  this  regime  cor- 
responds to  the  situation  In  which  the  viscous  sublayer  has  been  des- 
troyed. Thus  A^£j  is  effectively  zero  and  the  Inner  region  damping 
terms  D , D_,  and  D [as  given  in  Eqs.  (41)  through  (43)]  tend  toward 

1.0.  In  Finson' s model,  which  is  postulated  to  be  valid  for  three- 
dimensional  roughness  elements  less  than  the  boundary- layer  thickness, 
the  Individual  elements  are  considered,  assuming  attached  parallel  flow 
approaching  the  elements.  The  roughness  elements  each  provide  a distri- 
buted drag  for  the  mean  momentum  equation,  and  each  wake  provides  a 
source  for  the  generation  and  dissipation  of  turbulent  kinetic  energy. 
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Each  element  (and  all  elements  are  Identical  In  this  model)  has  a 
diameter  D^(y)  which  Is  a function  of  the  distance  from  the  wall  and  an 
average  center-to-center  spacing 


Figure  2.  as  a function  of  f for  the 
transitional  roughness  regime. 


The  form  drag  on  roughness  elements  represents  a sink  term  for  the 
streamwise  momentum  equation.  Referring  to  the  following  sketch,  the 
drag  between  (y  - 6y/2)  and  (y  + 6y/2)  on  a single  element  is 

-Js  p D^(y)  6y  (64) 
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y + 
y - 


_2 

There  are  roughness  elements  per  unit  surface  area  so  that  the 
momentum  sink  per  unit  differential  volume  is 

-1/2  p u^  CDj.  -y~  (65) 

A* 

r 

where  the  recommended  value  of  is  0. 5 with  a roughness  element 
spacing  of  Finson  (Ref,  33)  suggests  that  for  sand  grain  roughness, 

hemispherical  elements  spaced  about  two  base  diameters  apart  yield 
meaningful  results.  Equation  (65)  represents  the  additional  term  which 
appears  In  Eq.  (2),  the  streamwlse  momentum  equation,  and  Is  utilized 
only  when  fully  rough  turbulent  flow  Is  being  calculated.  The  velocity 
fluctuations  behind  an  element  are  taken  as 

u"  Qu  u (66) 

where  is  of  0 (0.1).  The  turbulent  kinetic  energy  generation  per 
unit  differential  volume  is 


o 3 °r(y) 

Q 7 U-*  — ^ — 

Xu  M u 2 

*r 


For  use  in  the  IKET  equation,  this  term  Is  Integrated  over  roughness 
height,  k^,  resulting  in 

“ I ^u  P (68) 

r 
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The  term  Is  then  added  to.Eq.  (37)  and  represents  the  turbulent 
kinetic  energy  generation  which  occurs  because  of  the  wakes  of  the 
roughness  elements.  A dissipation  term  could  also  be  Included,  but 
Finson  found  that  this  had  little  ef f ect‘.  Section  2.  5 contains  a dis- 
cussion of  the  interpretation  and  Importance  of  the  term  E^,  which 
appears  in  the  IKET  equation.  The  skin  friction,  or  more  properly,  the 
flow  resistance  Is  taken  as  the  sum  of  the  Integrated' form  drag  of  the 
elements  and  the  pure  frictional  force  applied  to  the  whole  surface;  - 
l.e. , ' ; ■■  ■ 


Furuya,  Miyata,  and  Fujlta  (Ref.  35)  express  the  data  from  their  experi- 
mental study  of  rough  wall  flows  In  such  a manner. 

All  of  the  roughness  examined  In  this  paper  was  assumed  to  be  of 
the  sand  grain  type.  This  is  obviously  an  optimistic  assumption  which 
greatly  facilitates  analysis  but  does  not  conform  to  the  roughness 
element  geometry  expected  In  aerodynamics.  For  other  than  sand  grain 
roughness  the  results  of  Dlrllng  (Ref.  36),  shown  in  Fig.  3,  are  used  to 
relate  various  roughness  geometries  to  an  effective  sand  grain  size,  and 
that  effective  sand  grain  size  Is  used  in  the  rough-wall  analysis. 

Reda,  Ketter,  and  Fan  (Ref.  37)  point  out  situations  In  which  the 
Dlrllng  correlation  should  be  used  with  caution. 

2^  IKET  EXTERNAL  SOURCE  TERMS 

The  IKET  equation  as  given  by  Eq.  (37)  contains  the  quantities  E 
and  E^,  which  represent  sources  of  turbulent  kinetic  energy  in  terms  of 
a global  energy  balance.  The  IKET  equation  In  effect  assimilates  the 
source  terms  E and  E^  by  varying  the  value  of  the  outer  region  length 
scale  1^,  namely  Y,  Local  values  of  quantities  such  as  the  turbulent 


31 


AEDC-TR-77-96 


2-D  Rod  Elements 


2-D  Wavy  Surface 


■ Bettermann 
4 Liu  et  ai. 

♦ Streeter 

3-D  Elements  - Schlichtinq 

o Spheres 
o Spherical  Segments 
^ Cones 

□ Cubes  - Koloseus 


• Streeter 
•"  Haughton 

3-D  Liquid  Patterns 
Cohen 


Recommended 
Correlation  Curve 

kg/kg  " 0.0164 

For  A < 4.93 

"139  A‘l-^ 

'For  A > 4.93 


Figure  3.  Effective  sandgrain  roughness  correlation. 

kinetic  energy  and  the  turbulent  shear  stress  are  then,  determined  by  the 
turbulence  model  using  the  distributions  of  the  length  scales  and 
and  the  structural  scale  a.j.  i 


Physically,  the  various  terms  which  comprise  the  quantity  E defined 
by  Eq.  (19)  represent  turbulence  source  terms  resulting  from  distur- 
bances imposed  upon  the  boundary  layer  by  the  free  stream.  The  quantity 
E is  actually  the  sum  of  two  contributions,  the  first,  , . . 
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ji/2  (p  u d6/dx  - p v)J  ^ (70) 

representing  the  free-stream  turbulence  entrained  by  the  boundary 
layer,  and  the  second, 

- 1/2^  - 1/2  p-  q^  (71) 

representing  the  direct  absorption  of  free-stream  acoustic  energy  by  the 
boundary  layer.  Both  constituents  of  E then  refer  to  quantities  at  the 
outer  edge  of  the  boundary  layer. 

The  quantity  E^,  defined  by  Eq.  (68),  represents  the  turbulent 
kinetic  energy  generation  related  to  the  wakes  of  individual  roughness 
elements  and  has  much  the  same  effect  on  the  IKET  equation  as  Edoes. 

But  unlike  E,  which  is  characterized  by  values  at  the  edge  of. the 
boundary  layer,  depicts  the  turbulent  kinetic  energy  generation  (by 
the  roughness  element  wakes)  Integrated  across  the  roughness  element 
geometry.  The  form  of  Eq.  (68)  Is  clearly  Indicative  of  the  Integrated 
nature  of  E^  as  well  as  of  the  dependence  upon  local  conditions  and  the 
geometric  shape.  Simply  stated,  free-stream  turbulence,  free-stream 
acoustic  disturbances,  and/or  roughness  element  wake  turbulence  serve  to 
provide  a small,  source  .term  for  the  IKET  equation. 

It  Is  this '‘source  which  mathematically  "triggers"  the  transition 
process  from  laminar 'to  turbulent  flow.  McDonald  and’ Fish  (Ref.  12)  and 
Shamroth  and  McDonald  (Ref.  13)  present  good  discussions  concerning 
initiation  of  transition  based  on  the  IKET  formulation.  However,  Eq'. 

(25)  for  the  turbulent  shear  stress  Is  strictly  applicable  only  under 
conditions  of  negligible  or  small  levels  of  free-stream  turbulence. 
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Hence,  the  present  IKET  formulation  is  restricted  to  negligible  or  small 
free-stream  turbulence  levels.  A recent  paper  by  McDonald  and 
Kreskovsky  (Ref,  14)  shows  how  this  restriction  may  be  removed. 

The  capability  to  assimilate  additional  energies  can  hardly  be 
overemphasized  as  it  allows  the  transition  process  to  be  computed  in  a 
continuous  fashion  from  completely  laminar  to  fully  developed  turbulent 
flow.  Within  the  constraints  of  the  current  analysis,  transition  can  be 
triggered  by  absorption  of  incident  acoustic  energy  and/or  roughness- 
element-generated  turbulent  kinetic  energy.  These  points  will  be  exam- 
ined further  in  Section  3.0. 

2.6  NUMERICAL  SOLUTION  OF  THE  COMPRESSIBLE  TURBULENT 
BOUNDARY-LAYER  EQUATIONS 

The  governing  compressible  time- aver aged  bound ary- layer  equations, 
either  two-dimensional  or  axisymmetrlc,  as  given  by  Eqs.  (1)  through 
(4),  as  well  as  the  IKET  equation  defined  by  Eq.  (37),  are  transformed 
to  Illlngworth-Levy  variables  and  formulated  into  linearized  finite- 
difference  equivalents;  the  finite-difference  forms  are  solved  using  an 
Iterative,  marching,  implicit  integration  technique  using  trldlagonal 
matrices.  Full  details  of  the  numerical  approach  may  be  found  in 
Appendixes  III  and  IV  of  the  reports  by  Adams  (Refs.  38  and  39).  The 
digital  computer  code  is  written  in  FORTRAN  IV  with  double-precision 
arithmetric  for  execution  on  an  IBM  370/165;  core  storage  is  180K  bytes 
using  the  FORTRAN  H OPT  «=  2 compiler.  As  with  other  parabolic  boundary- 
layer  codes,  outer-edge  conditions  must  be  specified  from  a separate 
source,  an  inviscid  analysis  or  an  experimental  pressure  distribution, 
for  example.  The  code  features  two-dimensional  or  axisymmetrlc  flow, 
sharp  or  blunt  leading  edges  (or  noses),  variable  grid  spacing  in  both 
streamwise  and  normal  directions,  and  arbitrary  wall  temperature  dis- 
tribution (or  adiabatic  wall) . 
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It 

coupled 

namely 


and 


Newton's  method  as  described  by  Conte  and  de  Boor  (Ref.  40)  Is  used  to 
solve  the  coupled  system.  The  outer-region  length  scale  “ Y appears 
as  a quadratic  function  In  the  IKET  equation  whose  solution  Is  numeri- 
cally determined  using  Newton's  method  for  root  extraction  as  given  by 
Householder  (Ref.  41).  A global  Iteration  procedure  Is  necessary  In 
solving  for: the  outer-region  length  scale  Y because  this  value  Is  linked 
with  the-iinner-reglon  scales  through  the  requirements  of  functional 
continuity; 

- - J.C-n  ■■■■'  ■ 

3.0  RESULTS  AND  DISCUSSION 
3.1  INCOMPRESSIBLE  FLAT  PLATE  FLOW 

I'.One  of  the '.'classical  experimental  Investigations  often  used  for 
comparison  with  analytical  theories  is  that  of  Klebanoff  (Ref,  42), 
-performed  ■some'- 20"years;?ago  in  the  National  Bureau  of , Standards  4-1/2-ft 
low-speed.'iwlnd  tunnel;  s An  artificially  thickened,  turbulent  boundary 
layer  was  allowed  to  develop  on  a smooth  flat  plate  12  ft  long.  The 
scheme  of  artificially  thickening  the  turbulent  boundary  layer  was 
achieved  by  covering  the  first  two  feet  of  the  plate  with  No.  16  floor- 
sanding,  paper.  Boundary- layer  probe  surveys  were  made  10. 5 ft  froni  the 
leading  edge,  at  which  point  the  boundary- layer  thickness  was  about  3 
In,  At  this  location  the  virtual  origin,  assiimlng  a smooth  surface,  was 


should  be  noted  that  In.  the  Inner  region  and  a.|  form  a 
nonlinear  set  of  algebraic  equations  at  a given  y- location, 


■■•■•[■">(-^¥1)] 


(72) 


(73) 
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14.2  ft.  The  corresponding  length  Reynolds  number  based  on  values  of  x 
measured  from  the  virtual  origin  was  4.2  x 10^  with  a free-stream 
velocity  of  50  ft/sec.  > 

Presented  in  Fig.  4 is  a comparison  of  the  mean  velocity  distri- 
bution as  measured  by  Klebanoff , and  the  IKET-computed  mean  velocity 
distribution.  The  agreement  is  excellent.  Also  illustrated  in  Fig.  4 
are  the  locations  of  the  inner  and  outer  regions,  as  well  as  the  near- 
wall region.  The  computed  turbulent  shear  stress  distribution  presented 
in  Fig.  3 agrees  as  well  with  Klebanoff s data  as  did  the  mean  profile. 
The  IKET  wall  shear  stress  value  (denoted  by  *)  is  in  good  agreement 
with  the  classical  Squire-Young  (Ref.  43)  correlation.  The  near-wall 
region  where  the  turbulent  shear  stress  is  damped  to  zero  is  not  shown 
in  Fig.  5. 

Figure  6 presents  a comparison  between  the  IKET-computed  turbulent 
kinetic  energy  distribution,  the  measured  turbulent  kinetic  energy 


—o—E)(perimental  Measurements  by  Klebanoff  (Ref.  42) 
X Present  IKET  Calculation 


Figure  4.  Mean-velocity  distribution  for  fuiiy  developed 
incompressible  turbulent  boundary  layer. 


36 


AEOC-TR-77-96 


■e—  Experimental  Measurements  by  Kletanoff  (Ref.  42) 
X Present  IKET  Calculation 

Wall  Shear  Stress  from  Present  IKET  Calculation 


distribution,  and  a turbulent  kinetic  energy  distribution  computed  by 
Pinson  et  al.  (Ref,  44).  using  a multi-equation  turbulence  model.  In  the 
outer  region  the  IKET  calculation  Is  In  good  agreement  with  the  data, 
but  In  the ■ inner  region  and  the  near-wall  region  the  experimental 
measurements  are  different  in  both  character  and  numerical  value  from 
the  IKET  calculation  and  the  more  sophisticated  approach,  of  Pinson  et 
al.  The  experimental  measurements  of  Klebanoff  as  presented. in  Fig.  6 
exhibit  every  small  amount  of  free-stream  turbulence  \*lch  is  reflected 
in  the  constant  turbulent  kinetic  energy  distribution  for  values  of  y 
greater  than  the  boundary- layer  thickness.  The  IKET  analysis  as  exam- 
ined in  this  report  assumes  zero  free-stream  turbulence  at  the  edge  of 
the  boundary  layer.  The  paper  by  McDonald  and  Kreskovsky  (Ref.  14) 
shows  how  the  restriction  of  zero  free-stream  turbulence  may  be  removed. 
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Figure  6.  Distribution  of  turbulent  kinetic  energy  in  fully 

developed  incompressible  turbulent  boundary  layer. 


3.2  COMPRESSIBLE  TURBULENT  BOUNDARY-LAYER  RELAMINARIZATION 


The  prediction  of  relamlnarlzatlon  phenomena  is  one  of  the  most 
stringent  tests  of  the  Innate  "physics"  of  turbulence  models.  Relami- 
narlzatlon  is  basically  a reversion  from  turbulent  to  laminar  boundary- 
layer  flow  caused  by  sustained  flow  acceleration  effects  such  as  occur 
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In  the  subsonic  , contraction  section  of  a nozzle.  One  of  the  best  docu- 
mented experimental  Investigations  of  compressible  relamlnarlzatlon  Is 
that  reported  by  Nash-Webber  (Ref.  45).  In  this  work  an  Instrumented 
flat  plate  was  opposed  by  a variety  of  upper-wall  profiles;  the  profiles 
were  chosen  to  Impose  various  pressure  gradients  on  the  flat-plate 
turbulent  boundary  layer.  Figure  7 Illustrates  the  edge  velocity  and 


Figure  7.  Velocity  atxl  acceleration  parameter  distribution 
for  Nash-Webber  Nozzle  A. 


the  free-stream  acceleration  parameter  K for  the  particular  upper  wall 

■ ■ ' iii  \ ‘ • 

(Nozzle  A)  used-^j.the  IKET-  calculation.  The  free-stream  acceleration 
parameter  Is  defined  as 


K = 


w 


Pw 


dU 


dx 


(74) 


where  the  w subscript  denotes  wall  conditions  and  the  e subscript  de- 
notes boundary- layer  edge  conditions.  The  Importance  of  this  parameter 
Is  Illustrated  by  Fig.  8,  taken  from  Ref.  45,  which  shows  that  the ' 
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Figure  8.  Turbulent-laminar  transition  boundary,  for 
adiabatic  wall  shear  layer. 


numerical  value  of  K can  be  used  as  an  Indicator  for  the  probable  oc- 
currence of  relamlnarlzatlon  provided  the  momentum  thickness  Reynolds 
number  based  on  edge  conditions  is  sufficiently  low.  Presented  in  Fig. 

9 are  calculated  results  from  the  present  IKET  analysis  and  the  experi- 
mental data  of  Nash-Webber.  Calculated  fully  laminar  and  fully  turbu- 
lent  results  are  also  shown  for  comparison  purposes.  The  present  IKET 
results  are  in  good  agreement  with  the  experimental  data.  The  local 
momentum  thickness  Reynolds  number  does  not  attain  a leunlnar  value  until 
X ^ 35  in.,  which  is  well  after  the  local  skin  friction  attains  its 
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laminar  value  (x  31  In.).  This  Is  due  to  profile  relaxation  effects 
of  the  boundary- layer  upstream  history,  a feature  not  possessed  by 
conventional  mixing-length  analyses.  Kreskovsky,  Shamroth,  and  McDonald 
(Ref.  46)  used  a similar  IKET  analysis  In  their  parametric  study  of 
relaminarlzatlon.  They  also  examined  the  Nash-Webber  Nozzle  A flow  and 
reported  good  agreement  between  their  IKET  analysis  and  experimental 
measurements. 


18  20"  22  24  26  28  30  32  34  36  38  40  42 
X.  in. 


Figure  9.  IKET  calculation  of  turbulent  boundary-layer 
, , . , .relaminarization. 


In  Fig.  9 the  equation 


(75) 
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simply  Indicates  the  IKET  extended  mixings length  expression,  which  Is 
equivalent  to  the  conventional  mixing- length  ratio  of  the  outer-region 
mixing  length  to  the  boundary- layer  thickness.  Conventional  mixing- 
length  analyses  nominally  assign  0.09  to  this  ratio,  but  the  IKET 
extended  mixing-length  approach  computes  the  ratio  In  the  streamwlse 
direction.  The  value  of  initially  characteristic  of 

turbulent  flow,  but  as  the  large  values  of  K are  approached  (x  as  28  in.) 
the  turbulent  boundary  layer  begins  to  undergo  relamlnarlzatlon,  which 
Is  reflected  in  the  decreasing  value  of  (A/5]  At  the  x- location 

where  the  skin  friction  approaches  the  laminar  value,  the  IKET  computed 
value  of  jKEi;  approaches  zero,  effectively  resulting  in  laminar 

flow.  Figure  9 also  points  out  that  the  value  of  starts  to 

decrease.  Indicating  the  beginning  of  relamlnarlzatlon,  well  before  any 
effect  on  the  skin  friction  is  felt.  This  Indicates  that  the  outer- 
region  length  scales  are  being  altered  by  the  pressure  gradient  well 
before  any  effects  are  felt  by  the  wall.  Thus  the  IKET  analysis 
possesses  something  akin  to  the  "precursor"  effect  (See  Section  3.4)  In 
relamlnarizlng  turbulent  boundary- layer  flows.  Figure  9 also  shows  that 
the  local  momentum  thickness  Reynolds  number  based  upon  boundary- layer 
edge  conditions  is  low  (about  2,000  at  the  x = 18-ln.  location);  there- 

‘‘‘f  f 

fore,  according  to  Fig.  8,  the  upstream  turbulent  boundary  layer  is 
Indeed  a candidate  for  relamlnarlzatlon,  provided  that  the  acceleration 
parameter  K Is  of  sufficient  magnitude.  By  comparing  the  results  of 
Fig.  9 with  the  acceleration  parameter  K given  in  Fig.  8,  one  can  see 
that  relamlnarlzatlon  occurred  over  the  region  where  the  numerical  value 
of  K was  greater  than  4 x 10  . The  crosshatched  region  in  Fig.  8 cor- 
responds to  the  current  numerical  example  of  the  Nash-Webber  Nozzle  A , 
and  confirms  the  basic  limits  of  validity  for  relamlnarlzatlon  as  given 
by  Kline  and  Launder  in  Fig.  8.  An  additional  example  of  relamlnarlzatlon- 
llke  behavior,  a rough  hemisphere-cylinder,  will  be  examined  in  Section 
3.7. 
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3.3  COMPRESSIBLE  TURBULENT  BOUNDARY  LAYERS  IN  ADVERSE 
AND  FAVORABLE  PRESSURE  GRADIENTS 


One  of  the  best  experimental  studies  concetmlng  compressible 
turbulent  boundary- layer  flows  under  nonequilibrium  conditions  is  that 
reported  by  Lewis*  Gran,  and  Kubota  (Ref.  47)  and  Gran,  Lewis,  and  Kubota 
(Ref.  48).  Details  of  the  AEDC  von  Karman  Gas  Dynamics  Facility  (VKF) 
Supersonic  Wind  Tunnel  (A)  experiment  may  be  found  in  the  report  of  Hahn 
and  Lutz  (Ref.  49).  The  wind  tunnel  model,  shown  in  Fig.  10,  consisted 
of  two  parts:  an  outer,  hollow  shell  whose  walls  can  be  cooled,  and  an 


Leading-E 
I nteractlon 
Shock 


p Compression 

\ ^Expansion  Fan 

\ \ r~Turbulent 


Boundary  Layer 


-Hollow 

Pressure 

Generating 

Centerbody 


Elnstr’umented  \ \ 


Figure  10.  Gran,  Lewis,  and  Kubota  wind  tunnel 
model  schematic. 

inner,  pressure-generating  body.  The  inner,  pressure-generating  body  was 
contoured  so  that  both  compression  and  expansion  regions  were  imposed  on 
the  outer  shell.  All  boundary- layer  measurements  (temperature,  pressure, 
heat  flux,  and  profile  surveys)  were  performed  on  the  Inner  surface  of 
the  outer  shell.  Nominal  free-stream  conditions  were 


M =3.98 

flO 

Re  /in,  = 5.18  x 10^ 
00 

T = STS^R 

0,“ 
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with  model  wall  temperature  ratios  of 

T /T  “ 0.9175  (hot  wall) 
w o 

T /T  “ 0.50  (cold  wall) 
wo,®. 

As  shown  in  Fig.  11,  the  centerbody  was  such  that  a zero  pressure 
gradient  initially  existed  (first  13.5  in.  of  the  model);  this  was 
followed  by  an  adverse  pressure  gradient  (from  13.5  to  18. 5 in  of  the 
model),  and  finally  by  a favorable  pressure  gradient  (over  the  remainder 
of  the  model).  The  severe  adverse  pressure  gradient  Increased  the 
surface  pressure  by  a factor  of  approximately  9 in  a streamwlse  dis- 
tance of  only  5 in.  Figure  11  also  Illustrates  the  Hach  number  varia- 
tion along  the  model  (analytical  expressions  are  given  in  Refs.  47  and 
48). 


— Present  I KET Calculation  Input 

Experimental  Data  from  Gran,  Lewis,  and  Kubota  (Ref.  48) 


inner  wall  of  Gran,  Lewis,  and  Kubota  iriodel. 


Comparison  of  the  IKET  calculations  and  the  experimental  data  are 
presented  in  Figs.  12  through  15  with  respect  to  displacement  thickness, 
momentum  thickness,  skin  friction,  and  Stanton  number  distributions. 
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- — Present  1KET  Calculation  . 


Figure  12.  Displacement  thickness  distribution  for  hot-wall 

condition  for  Gran,  Lewis,  and  Kubota  edge  conditions. 


X,  in. 


Figure  13. ' Momentum  thickness  distribution  for  hot-wall 

condition  for  Gran,  Lewis,  and  Kubota  edge  conditions. 
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O 


Experimental  Data  from  Ref.  42 
IKET  Calculations  Using  Lag  Model 


X,  in. 


Figure  14.  Skin  friction  distribution  for  hot-wall  condition 
for  Gran,  Lewis,  and  Kubota  edge  conditions. 


Present  IKET  Calculation 


Figure  15.  Stanton  number  distribution  for  cold-wall  condition 
using  Gran,  Lewis,  and  Kubota  edge  conditions. 
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respectively.  Except  for  the  momentuia  thickness  results,  the  agreement 
between  the  experiment  and  calculations  Is  reasonably  good.  As  can  be 
seen  from  Fig.  14,  the  variation  of  with  pressure  gradient  using  • 

the  lag  model  discussed  earlier  resulted  In  noticeably  better  IKET 
calculations  than  did  the  assumption  of  = 26.0  Irrespective  of  the 

pressure  gradient. 

Using  experimental  data  as  a basis,  Johnson  and  Kaufman  (Ref.  50) 
suggest  that  the  peak  heating  ratio,  (for  Interaction  of  shock  waves 
with  turbulent  boundary  layers),  can  be  estimated  In  terms  of  the  peak 
pressure  ratio,  as 


H 


pk 


'V 


(76) 


where  the  best  hypersonic  correlation  is  achieved  by  using  n = 0.85.  ''In 
Eq.  (76),  the  peak  heating  ratio,  is  defined  as  the  peak  heating 

rate  for  the  Interaction  flow  divided  by  the  heating  rate  for  undls- . 
turbed  flow  at  the  same  location,  and  the  peak  pressure  ratio, 
defined  as  the  peak  p/p^  for  the  Interaction  flow  divided  by  p/p^  for 
the  undisturbed  flow  at  the  same  location.  For  the  present  conditions, 
Ppj^  =9.1  (from  Fig.  11),  and  thus  = 6.53,  which  is  in  good  agree- 
ment with  the  results  of  Fig.  15. 


It  Is  germane  at  this  point  to  examine  the  variation  of  the  outer- 
region  length  scale. as  computed  using  the  IKET  analysis.  Figure  16 
presents  the  computed  hot-  and  cold-wall  boundary- layer  thicknesses  and 
outer-region  length  scales.  The  severe  adverse  and  favorable  pressure 
gradients  really  have  little  effect  on  the  outer-length  scale,  and  thus 
the  conventional  mixing-length  turbulent  boundary- layer  analyses  [based 
on  constant  outer-length  scale,  for  example  X/6  = 0.09  In  the  Fatankar- 
Spalding  (Ref.  10)  approach)  should  yield  relatively  accurate  results 
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for  this  condition.  Such  Is  indeed  the  case,  as  the  classical  analyses 
of  Reeves  (Ref.  51)  and  Cebecl  and  Smith  (Ref.  8)  report  good  agreement 

with  the  data  of  Refs.  47  and  48.  In  Fig.  16,  as  In  Fig.  9,  the  ex- 

' 

presslon  /C^/2a^  Y/6  represents  the  ratio  of  the  IKET  version  of  the 
conventional  mixing-length  outer-region  length  scale  to  the  boundary- 
layer  thickness. 

IKET  Calculation  for  Hot-Wall  Condition 

IKETCalculation  for  Cold-Wall  Condition 


0.08 
a45 
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0.35 
0.30 
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X,  in. 

Figure  16.  Outer  length  and  boundary-layer  thickness  distributions  ' ' 

for  Gran,  Lewis,  and  Kubota  edge  conditions.  iV-.i  ^ . 

The  highly  respected  klnetlc-energy-o£-turbulence  approach  by 
Bradshaw  and  Ferriss  (Ref.  21)  fails  dramatically  for  this  flow  situa- 
tion.  In  a later  paper  Bradshaw  (Ref.  52)  derived  an  empirical  bulk 

' ..if/ j a 

compression  or  dilatation  correction  formula  (essentially  an  "extra”  i. 
strain  rate)  which,  when  added  to  the  original  analysis,-  gave  good 
agreement  for  the  skin  friction  in  the  adverse  pressure  gradient  region. 
However,  this  modified  approach  still  failed  to  adequately  predict  the 
effects  of  the  favorable  pressure  gradient  region. 
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3.4  COMPRESSIBLE  BOUNDARY-LAYER  TRANSITION  ON  A 
SHARP  FLAT  PLATE 

The  comparisons  made  up  to  this  point  have  been  concerned  only  with 
the  prediction  of  turbulent  boundary  layers.  Boundary- layer  transition 
and  the  transitional  predictive  aspects  of  the  IKET  analysis  will  now  be 
examined.  The  classical  hypersonic  sharp  flat-plate  boundary-layer 
transition  experiment  of  Deem,  Erickson,  and  Murphy  (Ref.  53)  will  be 
used  as  a comparison.  Although  results  were  also  obtained  from  the 
AEDC-VKF  Hypersonic  Wind  Tunnel  (C)  and  the  AEDC-VKF  Supersonic  Wind 
Tunnel  (A),  only  the  Mach  No.  8 AEDC-VKF  Hypersonic  Wind  Tunnel  (B) 
conditions  will  be  considered.  The  free-stream  conditions  are 

M =8.09 

“ fi 

Re  /In,  = 0.25  x 10 
00 

T = 

0,“  ’ 

with  outer-edge  conditions 

M =7.40 

® f. 

Re  /in.  = 0.235  x 10 
e 

and  a model  wall  temperature  ratio  T /T  of  0.80.  Transition  was 

w o,« 

experimentally  determined  based  on  surface  pitot  pressure  measurements. 

The  experimental  studies  of  Pate  and  Schueler  (Ref.  54)  established 
that  boundary- layer  transition  on  models  in  supersonic  and  hypersonic 
wind  tunnels  is  dominated  by  radiated  aerodynamic  noise  from  the  tunnel 
sidewall  turbulent  boundary  layer.  The  recent  text  by  Goldstein  (Ref. 

55)  treats  on  a readable  analytical  basis  the  generation  and  propagation 
of  sound  waves  for  many  aerodynamic  circumstances.  The  radiated  aero- 
dynamic noise  from  the  supersonic  wind  tunnel  sidewall  turbulent  bound- 
ary layers  Is  normally  characterized  by  the  level  of  the  root-mean- 
square  (RMS)  radiated  pressure  fluctuations.  Recent  measurements  by 
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Lademian  (Ref.  56)  Indicate  that  for  the  aforementioned  free-stream 
conditions  the  RMS  radiated  pressure  fluctuations  in  the  AEDC-VKF  Tunnel 
B are  approximately  3.5  percent  of  the  free-stream  static  pressure. 

Thus  the  E term  In  the  IKET  equation  [Eq.  (37)]  represents  the  incident 
acoustic  energy  (In  the  form  of  radiated  pressure  fluctuations)  absorbed 
by  the  boundary  layer.  It  Is  this  source  term  which  triggers  transition. 

The  very  difficult  problem  of  receptivity  as  defined  by  Morkovin 
and  discussed  by  Reshotko  (Ref.  57)  is  Intimately  associated  with  the 
theoretical  determination  of  E.  In  this  analysis  E Is  evaluated  via 
experimental  data,  and  no  effort  is  made  to  analytically  determine  the 
absorbed  incident  acoustic  energy  by  considering  the  incident  spectrum 
and  its  signature  within  the  boundary  layer.  This  is  accomplished  by 
selecting  the  value  of  E which  places  the  computed  transition  point  at 
the  desired  location  using  the  appropriate  criterion,  In  this  example 
the  maximum  value  of  the  surface  pitot  measurement.  Presented  In  Figs. 
17,  18,  and  19  are  numerical  results  from  the  IKET  analysis  which  were 
Initiated  via  a laminar  similar  solution  at  the  sharp  leading  edge  and 
Integrated  downstream  using  a constant  value  for  the  E source  term.  The 


X,  In, 

Figure  17.  Laminar,  transitional,  and  turbulent  displacement 
thickness  distribution. 
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Figure  18.  Laminar,  transitional,  and  turbulent  momentum 
thickness  distribution. 

experimental  and  computed  boundary- layer  displacement  thicknesses  are  in 
good  agreement.  The  transition  phenomena  are  more  clearly  seen  in  Fig, 

19,  where  the  calculated  outer-region  length  scale  starts  at  zero  at  the 
leading  edge  and  Increases  downstream,  reflecting  the  transition  process. 
At  the  X = 11-in.  station,  which  is  the  location  of  minimum  skin  friction, 

fl 

the  outer-length  scale  ratio  /C  /2a  y/6  has  attained  a value  of  0.034, 

which  is  approximately  40  percent  of  the  maximum  value  of  0. 08,'  attained 
when  transition  to  turbulence  is  complete.  Hence,  the  IKET  analysis 
predicts  substantial  values  of  the  outer-region  length  scales  before  the 
surface  feels  any  effects  of  transition.  Further  discussion  of  the 
length  scale  behavior  during  transition  will  be  found  later  in  this 
section.. 

Also  shown  in  Fig.  .(19  are  IKET- calculated  results  for  cold-wall 
conditions.  The  influence  of  the  cold  wall  is  to  slightly  displace 
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downstream  the  minimum  skin-friction  location  and  the  outer  length 
scale.  Hence,  the  IKET  analysis  predicts  the  slight  stabilizing  effect 
of  a cold  wall  on  boundary- layer  transition.  This  agrees  with  the 
meager  experimental  results  of  Deem  et  al.  (Ref.  53)  and  with  the  gen- 
erally accepted  transition  trends  [whereas  wall  cooling  stabilizes  the 
laminar  boundary  layer  relative  to  transition  (Ref.  57)1.  "Hie  cold-wall 
Influence  on  turbulent  skin  friction  (l.e. , an  Increase  over  the  cor- 
responding hot-wall  value)  should  also  be  noted. 


Sharp  Flat  Plate,  AEDC-VKF  Tunnel  B 
Mg  *7. 40,  Reg/in.  *2.35x105 

I KET  Calculation  for  oo  “ 

I KET  Calculation  for  Ty/Tg^  jg  * 0. 40 


Figure  19.  Laminar,  transitional,  and  turbulent  outer  ^ , 

length  scale  and  local  skin-friction  distributions. 

* .'  ■ '■}  f' t 

Any  analytical  model  of  the  boundary- layer  transition  process,  under 
hypersonic  conditions  must  include  the  precursor  effect  as  discussed-  by 
Bushnell  and  Alston  (Ref.  58).  This  effect  is  characterized  by  the 
existence  of  appreciable  length  scales  (l.e.,  large-scale  disturbances) 
in  the  outer  region  of  compressible  boundary  layers  well  upstream  of  the 
nominally  accepted  transition  location.  Figure  20,  taken  from  Ref.  58, 
schematically  illustrates  the  phenomenology  of  the  precursor  effect. 
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Location -x^'  Indicates  the  point  of  the  Initial  breakdovm  of  laminar, 
flow.  At  this  location  laminar  profiles  and  characteristics  prevail. 
Location  represents  the  point  at  which  turbulent  bursts  are  initiated 
In  the  outer  region  of  the  still  laminar  boundary  layer  and  the  point  at 
which  appreciable  growth  of  the  outer-region  length  scales  begins.  The 
location  at  which  transition  is  indicated  by  surface  measurables  such  as 
skin  friction  and  heat  transfer  is  represented  by  x^.  Since  turbulent 
bursts  are  observed  far  upstream  of  the  conventional  transition  location 


I nitlal  Breakdown  of  I nltiatlon  of  Rapid  I ncrease  with  x of 

Laminar  Flow  (In  Turbulence  Wall  Transport  (4,  Cf) 

Outer  Region)  "Bursts"  "Conventional"  Transition 

Location,  Profile  Already 
Distorted  in  Outer  Region 

Figure  20.  Schematic  of  the  precursor  effect  in  hypersonic 
boundary-laver  transition. 

' .*  >■  " 

where  the  wall  heat  transfer  or  skin  friction  deviates  from  the  laminar 
values,  transition  at  hypersonic  speeds  must  occur  initially  away  from 
the  wall  as  per , the  ’wrk  of  Potter  and  Whitfield  (Ref.  59)  and  Stainback 

'iv;;-,**;;  nr;-*?:;"’'  • i : ■ , ■ 

(Ref.  60).  Thus,  the  outer  portion  of  the  mean  flow  profiles  may 
already  be> distorted' by  turbulence  effects  at  the  conventional  transi- 
tion location.  This  was  shown  by  Fischer  and  Weinstein  (Ref.  61)  and 
is  contrary  to  the  usual  numerical  approach,  e.g.,  Adams  (Ref.-  62), 
Harris  (Ref.  63),  and  Kuhn  (Ref.  64),  where  the  intermittency  distri- 
bution is  anchored  to  the  conventional  transition  location  with  a fully 
laminar  profile.  This' precursor  effect  is  predicted  properly  by  the 
IKET  analysis  as  shown  in  Fig.  21,  which  gives  the  turbulent  shear 
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stress  profiles  for  three  streamwise  locations  along  the  plate,  namely 
6,  12,  and  18  in.  The  IKET  analysis  predicts  at  the  6- in.  location, 
which  is  far  upstream  of  the  surface- indicated  transition  location,  a 
sharp  "spike"  in  the  turbulent  shear  stress  near  the  outer  edge  of  the 
boundary  layer.  The  transitional  turbulent  shear  stress  profile  at  x ° 
12  in.  approaches  the  fully  turbulent  value  at  the  x ° 24- in.  location, 
especially  for  the  hot-wall  condition.  Thus,  in  the  outer  region  of  the 
boundary  layer  at  x = 12  in.  the  flow  is  turbulent,  and  hence  the  outer 
portions  of  the  boundary- layer  mean  flow  profiles  are  Indeed  distorted 
by  upstream  turbulence  effects. 

Boundary-Layer  Thickness  6,  in. 

X,  in.  Tvi)/To^  (j)  • 0. 40  Tvy/To^  o)  • 0. 80 


1/2  PeUe^ 


Figure  21.  Turbulent  shear  stress  profiles  as  an  indicator 
of  the  precursor  effect. 

Defining  the  distance  y^^  shown  in  Pig.  20  as  the  y location  \diere 
the  turbulent  shear  stress  attains  its  maximum  value  allows  the  results 
of  Fig.  21  for  the  x = 6- in.  location  to  be  tabulated  as 
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T /T 

y /6 

W O," 

■^cr 

0.40 

0.89 

0,80 

0.92 

which  Is  In  good  agreement  with  the  results  of  Potter  and  Whitfield 
(Ref.  59)  and  Stalnback  (Ref.  60).  Since  the  value  of  the  maximum 
turbulent  shear  stress  at  x ° 6 In.  Is  approximately  one  percent  - of  the 
maximum  turbulent  shear  stress  at  x =■  12  in..  It  Is  reasonable  to  pos- 
tulate that  the  IKET-calculated  concentration  of  turbulent  shear  stress 
near  the  edge  of  the  boundary  layer  represents  In  a time-mean  sense  the 
turbulent  "burst"  observed  experimentally  by  Potter  and  Whitfield  (Ref. 
59)  as  well  as  by  Fischer  and  Weinstein  (Ref.  61). 

The  edge  Mach  number  for  this  calculation  was  7.40,  a low  hypersonic 
Mach  number  which  exceeds  the  formal  limit  of  validity  of  Morkovln's 
hypothesis  (Ref.  27).  Nevertheless,  the  results  of  the  computations  for 
this  example  appear  satisfactory.  Indicating  that,  at  the  formal  limit 
of  validity,  "breakdown"  is  slow  rather  than  catastrophic.  The  Morkovln 
hypothesis  as  used  In  this  analysis  permits  structural  constants  as 
deduced  from  Incompressible  turbulent  flow  to  be  used  for  edge  Mach 
numbers  up  to  about  five.  Thus,  for  the  edge  Mach  number  (7.40)  used  In 
this  example,  spine  .turbulence  model  constants  may  In  reality  deviate 

■'L'.i  ' . , 

slightly  from  the^ assumed. values,  but  any  deviation  must  be  very  gradual 
at  this  low  hypersonic  Mach  number,  and  the  resulting  Impact  upon  the 
calculations  must  be  correspondingly  small. 

Figure  19  also  indicates  that  the  outer  region  length  scale  ratio 

y2  ■ ^ 

expressed  as  /C^/2a^  Y/&  possesses  a value  near  0.08.  The  usual  as- 
sumption made  in  mixing-length  hypotheses  is  that  the  outer-region 
length  scale  ratio,' X/ 5,  has  a value  of  0.09.  Previous  examples  con- 
sidered in  this  report. have  tended  to  confirm  the  vaiue’ of  0.  09.'  This; 
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example  raises  the  possibility  that  for  hypersonic  flows  the  outer- 
region  length  scale  ratio  may  decrease  from  the  usually  accepted  value 
of  0.09. 


The  determination  of  the  incident  acoustic  energy,  E,  absorbed  by 
the  boundary  layer  is  greatly  facilitated  by  the  use  of  IKET-generated 
plots,  of  which  Fig.  22  Is  Illustrative.  This  figure  presents  the 
relationship  between  Re  , the  edge  Reynolds  number  at  transition,  and 

^ 6yXC  ^ 

E , the  incident  absorbed  acoustic  energy  made  dimensionless  by  p U*^. 

6 . e e 

For  the  sharp  flat  plate  at  a given  Mach  number  and  wall  temperature, 
these  parameters  collapse  all  free-stream  Reynolds  numbers  and  values  of 
E onto  a single  curve  which  is  a straight  line  when  displayed  on  a log- 
log  plot.  Thus,  the  absorbed  incident  acoustic  energy  required  for 
transition  at  a given  location  can  be  readily  ascertained.  Implicit  in 
this  procedure  for  obtaining  the  required  incident  acoustic  energy 
absorbed  by  the  boundary  layer,  E,  is  that  the  receptivity  be  handled  in 
an  "effective"  manner.  Linear  compressible  stability  theory  as  devel- 
oped by  Mack  (Ref . 65)  can  be  used  to  examine  qualitatively  the  problems 
inherent  in  the  theoretical  calculations  of  E.  The  results  of  Mack  for 
a flat  plate,  as  reported  by  Reshotko  (Ref.  57),  show  that  in  the  Mach 
number  range  from  2.5  to  7 all  sound  frequencies  grow  rapidly  monotoni- 
cally  from  the  leading  edge,  thus  indicating  that  the  acoustic  energy 
absorbed  by  the  transitioning  boundary  layer  is  a function  of  the  RMS  '• 
radiated  pressure  fluctuations.  Above  a Mach  number  of  about  seven  the' 
laminar  boundary  layer  becomes  much  more  selective,  and  only  distur- 
bances of  discrete  frequencies  are  amplified.  This  leads  to  behavior 
in  which  the  energy  absorbed  by  the  boundary  layer  is  dependent  upon  " ■ 
both  the  disturbance  environment  spectral  energy  distribution  and  the 
detailed  structure  of  the  boundary  layer.  The  value  of  E as  determined 
from  graphs,  of  which  Fig.  22  is  typical.  Implicitly  contains  all  the 
nuances  of  receptivity  for  a single  given  case.  Nothing  additional 
about  E in  general  can  be  inferred. 
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Figure  22.  Transition  Reynolds  number  versus  absorbed 

acoustic  energy  for  a sharp-leading-edge  flat  plate. 

3.5  COMPRESSIBLE  BOUNDARY-LAYER  TRANSITION  ON  A 
SHARP  CONE 

. The  f and  results  for  axlsymmetrlc  boundary  layer 

on  bodies  of  revolution  will  be  for  sharp  cones.  The  same  turbulence 
njodel;7iiiay  be  u8ed,in,lK£T  analysis  of  axlsymmetrlc  boundary  layers  as 
was, used  in  the  I^T., analysis  of  planar  boundary  layers.  Presented  In 
Fig.; 23  tare  sharp-cone. transition  Reynolds  numbers  (based  on  edge 
conditions|)  versus. dimensionless  Incident  absorbed  acoustic  energy,  E^. 
The  resultStjShown-.  in  Fig,.  23  for  sharp  cones  are  very  reminiscent  of  the 
results  portrayed  in  Fig..  22  for  sharp  flat  plates. 


. The  results,- of  Pate  (Ref.  66)  for  sharp,  slender  cones  show,  that  ■ 
the  ratio. of  the . cone,  transition  Reynolds  number  to  the.  planar  transi- 
tion Reynolds  number  varies  from  about  2.5^at  a Mach  number  of,  3.0  to  , 
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Figure  23.  Transition  Reynolds  number  versus  absorbed 
acoustic  energy  for  a sharp  cone. 


approximately  1.0  at  a Mach  number  of  8.0.  These  results  (taken  from 
Fig.  11  and  Table  3 of  Ref.  66)  are  reproduced  in  Fig.  24,  where  the 
legend  delineates  the  meaning  of  the  open,  solid,  and  flagged  symbols. 
The  solid  line  in  Fig.  24  represents  the  IKET-predlcted  transition 
Re3rnolds  number  ratios  as  a function  of  the  edge  Mach  number.  The  solid 
line  was  generated  using  identical  edge  conditions  for  the  sharp  cone 
and  the  flat  plate  and  identical  values  of  the  absorbed  incident  acous- 
tic energy,  E.  The  results  of  Fig.  24  suggest  that  the  acoustic  energy 
absorbed  by  both  a sharp-cone  and  a flat-plate  boundary  layer  (both  with 
the  same  edge  conditions)  must  be  the  same. 

Deviations  begin  to  appear  at  a Mach  number  of  about  six,  which  is 
near  the  limit  of  validity  of  Jforkovin’s  hypothesis  (Ref.  27).  As 
discussed  in  Section  3.4,  the  formal  range  of  validity  of  some  turbu- 
lence model  constants  has  probably  been  exceeded  by  Mach  number  six,  but 
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the  rather  drastic  departure  of  the  IKET-predicted  results  for  the 
higher  Mach  numbers  Is  not  Indicative  of  what  must  be  a rather  gradual 
breakdown  In  the  turbulence  model.  Moreover,  previous  calculations 
reported  herein  have  Indicated  reasonable  results  for  low  hypersonic 
Mach  numbers.  Linear  compressible  stability  theory  offers  a possible 
explanation  for  the  anomaly  observed  in  Fig.  24.  As  pointed  out  In 
Section  3.4,  linear  compressible  stability  theory  Indicates  that  above  a 
Mach  number  of  about  7 (for  flat  plates  at  least)  the  laminar  boundary 
layer  becomes  selective  In  the  frequency  of  disturbances  amplified.  The 
IKET  analysis  as  applied  In  this  example  has  essentially  assumed  energy 
absorption  on  an  RMS  basis  (l.e.,  for  given  edge  conditions,  E values 
for  both  planar  and  conical  flows  were  taken  to  be  the  same).  However, 


Figure  24.  Cone  and  planar  transition  Reynolds  number 
ratios  as  a function  of  local  Mach  number. 


59 


AEDC-TR-77-96 


in  regions  where  the  growth  of  disturbances  is  highly  dependent  upon 
disturbance  frequency  and  detailed  boundary-layer  structure,  a broadband 
representation  such  as  this  should  not  be  successful.  Differing  recep- 
tivities of  planar  and  conical  boundary  layers  at  the  higher  Mach 
numbers  then  provide  a reasonable  explanation  of  the  anomaly  observed  at 
higher  Mach  numbers  in  Fig.  24.  An  Interesting  discussion  of  the 
relationship  between  linear  stability  theory  on  a cone  and  a flat  plate 
is  given  by  Mack  (Ref . 65) . 

Previous  results  discussed  have  established  the  relationship 
between  the  transition  location  and  the  absorbed  incident  acoustic 
energy,  but  no  mention  has  been  made  of  the  effects  (if  any)  of  absorbed 
incident  acoustic  energy  upon  turbulent  boundary  layers.  The  results 
presented  in  Fig.  25  show  the  effect  of  various  energy  absorptions  upon 
the  IKET-computed  skin  friction  for  typical  wind  tunnel  flow  conditions 

Ma)-3.0 


Figure  25.  Effect  of  acoustic  energy  absorption  on  a 
turbulent  boundary  layer. 
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over  a 5-deg  half-angle  sharp  cone.  The  solid  line  in  the  figure  was  ■ 
calculated  using  the  usual  IKET  transition  analysis,  in  which  an  acous- 
tic energy  absorption  is  specified  and  the  resulting  transition  process 
is  calculated;  the  dashed  line  was  computed  by  forcing  transition  to 
occur  instantaneously,  with  no  energy  absorption;  the  symbols  were  i ‘ ' 

computed  using  forced  Instantaneous  transition  and  values  of  energy'  < 
absorption  10  and  100  times  that  required  for  the  transition  calcula-  ■ 
tlon.  A decade  Increase  in  the  energy  absorption  has  almost  no  dls-  ■ 
cernible  effect  upon  the  calculated  skin  friction,  whereas  the  large 
multiple  of  100  gives  about  10  percent  increase  above  the  other  cal- 
culations. This  latter  amount  represents  almost  total  absorption  of  all 
acoustic  energy  incident  upon  the  boundary  layer  and  hence  can  bie  con-  - 
sldered  an  unrealistic  number.  Thus,  for  typical  calculations  the 
effect  of  , absorbed  acoustic  energy  upon  turbulent  boundary  layers  is 
small.  It  ;,can  therefore  be  inferred  that  in  typical  supersonic  wind 
tunnels  the  acoustic  environment  has  considerable  effect  upon  the  transi- 
tion location  but  little  effect  upon  the  turbulence  level  in  a fully  ' ' 

turbulent  boundary  layer, 

3.6  SURFACE  ROUGHNESS  IN  PLANAR  COMPRESSIBLE 
TURBULENT  BOUNDARY  LAYERS 

Roughness  effects  were  investigated  using  as  a basis  of  comparison 
the  results  of  Reda,  Ketter,  and  Fan  (Ref.  37),  which  were  taken  on  the 
flat  nozzle  wall , of -a  ..supersonic  nozzle  in  which  the  wall  opposite  was 
contoured.  A wide;  variety  of  .roughness  conditions  was  studied  over  a 
wide  range  of  free-stream  Reynolds  numbers.  Skin-friction  and  profile 
data  were  taken,  with  the-  skin-friction  data  being  obtained  from  a 
floating  balance.  Examined  in  this  report  are  the  data  for  smooth 
walls  and  sand  grain  roughness  of  80,  50,  and  24  grit.  This  range  of 
roughness  allows  all  three  regimes  of  rougl^wall  turbulent  boundary- 
layer  flow  to  be  examined. 
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Figure  26  contains  a summary  of  the  predictive  ability  of  the 
rough-wall  IKET  analysis.  The  smooth-wall  computed  results  and  experi- 
mental data  show  the  same  trend  with  increasing  Reynolds  numbers,  but 
the  predicted  level  is  slightly  low.  This  is  difficult  to  explain  but 
not  too  disturbing,  as  changes  in  free-stream  conditions  and/or  the 
location  of  the  boundary  layer's  virtual  origin  in  the  computed  results 
or  sensitivity  of  the  balance  in  the  experiment  could  have  caused  such 
an  error.  The  80-grit  data  are  quite  Interesting  in  that  initially  the 
roughness  Reynolds  number  was  very  close  to  the  "smooth"  limit;  hence. 


Flat-Plate  Turbulent  Boundary  Layer, 
.Mqo-2.9,  Adiabatic  Wall 


Figure  26.  Skin  friction  versus  for  ' 
rough-wall  flow. 


the  resulting  nearly  smooth-wall  value  of  skin  friction.  As  the  Rey- 
nolds number  Increased,  the  roughness  became  more  pronounced  with 
respect  to  the  viscous  sublayer,  causing  a thinning  which  is  reflected 
in  the  higher  value  of  skin  friction.  The  80-grlt  roughness  remained 
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continuously  In  the  transitional  roughness  regime,  thus  allowing  the 
Healzer  et  al.  (Ref  . 30)  approach  (see  Section  2.4)  tp,  be  used.  , The  50- 
grit  roughness  was  more  difficult,  as  at  the  low  free-stream  Reynolds  , 
numbers  the  roughness  was  transitional,  and  at  higher  free-stream  ..  . . 
Reynolds  numbers  the  roughness  was  fully  rough.  The  fully, rough  regime  . 
was  treated  using  the  approach  of  Finson  (Ref.  33.,  see  Section  2.4). 

The  fully  rough  calculations  used  the  suggestion  of  Finson  that  sand 
grain  results  could  be  adequately  predicted,  assuming  hemispherical 
roughness  elements  of  a height  c^racteristic  of,  the'  grit  specification 
(0.00225  ft  for  24-grlt,  for  example)  and  a center-to-center  spacing  of 
two  element  diameters.  The  results  obtained  were  considered  satisfac- 
tory, and  no  computer  experiment  was  undertaken  to  examine  various 
' configurations  of  roughness  element  shape  and  spacing.  The  24-grit  sand 
grain  roughness  was  large  enough  so  that  the  viscous  sublayer  was  com- 
pletely destroyed  at  all  the  free-stream  Reynolds  nxjmbers  of  Interest, 
and  the  flow  was  fully  rough,  thus  allowing  the  Finson  approach  to  be 
used  for  all  Reynolds  numbers. 

For  some  free-stream  conditions  the  roughness  was  such  that  the 
flow  was  initially  fully  rough,  but  as  the  viscous  sublayer  thickened, 
the  flow  reverted  into  the  transitional  regime.  No  completely  satis- 
factory method  of  .-interfacing  the  fully  rough/ transitional  regimes  was 
evolved,  but  results  appeared  reasonable  if  the  Finson  fully  rough 
analysis,  once  Invoked,  was  used  until  a local  roughness  Reynolds 
number  below  40.0  was  encountered.  The  Healzer  et  al.  (Ref.  30) 
transitional  roughness  analysis  was  then  applied. 

3.7  ROUGH-WALL  HEMISPHERE-CYLINDER  IN  HYPERSONIC  FLOW 

Hemisphere-cylinders  have  traditionally  been  used  as  standardized 
test  shapes  in  supersonic  and  hypersonic  aerodynamics  because  they  offer 

' ■:q  , • . , - - • . 

a wide  range  of  edge  pressures,  edge  temperatures,  and  pressure  gradients 
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Recently  the  AEDC-VKF  Hypervelocity  Wind  Tunnel  (F),  which  usea  nitrogen 
as  a working  medium,  was  used  to  investigate  the  effectiveness  of  sand 
blasting  model  surfaces  in  order  to  Induce  turbulent  flow.  A hemisphere- 
cylinder  model  with  a nose  radius  of  2.5  in.  was  Instrumented  with  heat- 
transfer  gages  and  sand  blasted  over  one  sextlle  as  shown  in  Fig.  27. 
Typical  run  conditions  in  Tunnel  F were  as  follows: 

M - 9.0 
00 

Re  /ft  - 3.40  X 10^  to  5.85  x 10^ 

00 

T - 19rR 
00 

T - 540“R 
w 

Nitrogen 


Figure  27.  Tunnel  F rough-wall  hemisphere-cylinder  model. 


Figure  28  presents  experimental  heat-transfer  data  for  the  test  condi- 
tions delineated  above.  The  solid  symbols  represent  data  taken  on  the 
sand-blasted  sextlle  of  the  model,  whereas  the  open  symbols  represent 
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data  taken  on  the  "smooth"  portion  of  the  model.  For  reference,  theo- 
retical computations  made  for  laminar  flow  and  for  fully  turbulent  flow 
using  a conventional  mixing-length  approach  (BLIMP,  Ref.  67)  are  shown. 


Figure  28.  Rough-wali  hemisphere-cylinder  in 
hypersonic  flow. 


The  experimental  data  Initially  Indicate  heat-transfer  rates  in  excess 
of  the  fully  turbulent  calculations.  The  discrepancy  is  caused  by  the 
enhanced  heat-transfer  rate  to  the  rough  surface,  an  effect  not  accounted 
for  in  the  fully  turbulent  BLIMP  calculations.  The  real  region  of 
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interest  occurs,  however,  between  the  x/R^  = 1.0  and  x/R^^  = 1,5  stations. 
In  this  region  the  measured  heat-transfer  rates  drop  below  the  predicted 
fully  turbulent  values  and  approach  the  laminar  predicted  values.  Thus 
the  data  indicate  that  relamlnarizatlon,  or  something  closely  akin,  is 
occuring  in  this  region.  Figure  29  shows  the  edge  pressure  ratio  (p/p^), 
the  edge  Mach  number  (M^),  and  the  acceleration  parameter  K [defined  by 
Eq.  {74)1  as  computed  using  the  AEDC  version  of  the  NASA  Ames  character- 
istics computer  code  (Ref.  68)  with  the  nose  solution  generated  by  the 
Aungier  program  (Ref.  69).  Near  the  =1.5  location,  K Increases 
dramatically  to  a value  near  1 x 10~^.  The  local  maximum  value  of  K at 
this  location  indicates  that  according  to  the  criterion  of  Fig.  8, 
relaminarlzation  is  just  possible  irrespective  of  Re  Actually  the  K 
value  falls  so  close  to  the  threshold  boundary  indicated  in  Fig,  8 that 
it  is  difficult  to  assess  its  meaning;  1. e, , relamlnarizatlon  may  or  may 
not  take  place,  in  spite  of  the  high  values  of  K near  the  x/R^j  =0.5 
region,  relamlnarizatlon  did  not  occur,  as  roughness  effects  were 
dominant.  Calculations  were  made  for  this  model  in  which  a smooth-wall 
flow  was  tripped  turbulent  initially  and  then  computed  downstream  from 
that  location  using  the  IKET  analysis.  Relamlnarizatlon  with  the  atten- 
/"2 

dant  approach  of  /C^/2a^  Y/5  to  near-zero  values  was  always  observed 
prior  to  x/R^j  =>  0.25.  Therefore,  the  presence  of  roughness  can  alter 
the  threshold  boundary  for  relaminarlzation  as  shown  in  Fig.  8.  Of 
course  relaminarlzation  could  still  occur  at  x/Rjj  near  1.50  since  k^/6 
decreases  from  0.2  near  the  x/K^  =0.5  location  to  0.02  near  the  x/R^^  = 
1.50  location,  indicating  that  the  relative  roughness  hAs  decreased  by 
an  order  of  magnitude.  ' 

• :a  T ion  v v ' S; . 

Figure  28  presents  the  results  of  the  IKET  rougl^wall  analysis  for 
the  conditions  of  Fig.  29.  The  IKET  analysis  with  E,  the  incident 
acoustic  energy  absorbed  by  the  boundary  layer,  taken  as  zero  effec- 
tively and  realistically  computes  both  the  roughness- Induced  heat- 
transfer  augumentatioh  in  the  x/R^^  = 0.5  region  and  the  relamlnarlza- 
tlon-like  process  near  the  x/R^^  = 1.5  region.  The  least  pleasing  result 
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is  observed  near  the  = 1.80  station,  where  the  predicted  heat- 

transfer  ratio  is  somewhat  above  the  measured  values.  No  explanation  is 
apparent,  as  the  remainder  of  the  cylindrical  portion  Is  adequately 
predicted.  Fig.  30  reveals  a moat  intriguing  result,  namely  that  the 


Me 


K 


x/Rn 


Figure  29.  Edge  conditions  and  acceleration  parameter  for 
M.  = 9 heroisphere«ylinder. 

outer-region  length  scale  /C^/2a^  Y/6  exhibits  turbulent,  not  laminar 
behavior.  Thus,  relaminarlzation  in  the  sense  of  the  example  in  Section 
3,2  did  not  occur  near  the  x/R^  “ 1*5  region  of  the  hypersonic  rough- 
wall  hemisphere-cylinder.  Additional  computations  were  made  in  which 
the  baseline  turbulence  model  values  of  the  parameters  k and  were 
assumed.  These  computations  showed  no  lamlnar-llke  behavior  in  the 
heat-transfer  rate  about  the  = 1.5  region.  Hence,  the  IKET  anal- 
ysis indicates  that  the  laminar-like  behavior  observed  near  x/R^^  =»  1,50 
in  the  rough-wall  hemisphere-cylinder  data  from  Tunnel  P is  caused  not 
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by  a reversion  to  laminar  flov  but  by  the  acceleration  (large  favorable 
pressure  gradient)- Induced  growth  of  the  viscous  sublayer.  The  value  of 
A^ff,  which  Is  a measure  of  the  viscous  sublayer  thickness  (see  Section 
2.4),  Is  computed  to  be  as  large  as  40  In  the  region  of  Interest.  This 
represents  a very  large  Increase  In  the  effective  thickness  of  the 
viscous  sublayer,  as  was  small  near  the  fully  rough/ transitional 

roughness  Interface.  The  Interpretation  of  Kays  and  Moffat  (Ref.  29), 
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in  which  relaminarization  is  attributed  to  the  viscous  sublayer's  . 
overwhelming  the  remainder  of  the  boundary  layer  (very  large  • j, . 

values  of  » thus  suggests  that  the  favorable  pressure  gradient  in  : 

the  proximity  of  = 1.50  was  not  strong  enough  to  result  in  a vis-  ■ 

cous  sublayer  thickness  approaching  the  boundary- layer  thickness. 
Because  the  viscous  sublayer  thickness  remained  relatively  thin  * 

40) , the  outer-region  length  scale  as  computed  using  the  IKET  analysis 
did  not  show  any  tendency  to  approach  zero.  Hence  at  the  =1.50 

region  the  flow  could  be  viewed  as  verging  on  relaminarization.  The  K 
value  of  near  1 x 10  ^ and  the  IKET-computed  value  of  Re  of  700  at 

e j y 

x/Rjj  = 1.5  fall  so  close  to  the  threshold  curve  as  given  in  Fig.  8 that 
the  conclusion  of  impending  relaminarization  is  reinforced. 


The  outer-region  length  scale  ratio  as  shown  in  Fig.  30  never 
exceeds  a value  of  0.086  and  appears  to  be  approaching  a value  of  0.080. 
As  in  previous  low  hypersonic  examples  (Section  3.4)  the  outer-region 
length  scale  ratio  is  less  than  the  usually  accepted  0.09.  The  work  of 
Bushnelly-  Cary,  and. Holley  (Ref.  70)  established  that  the  value  of  the 
outer-region  length  scale  ratio,  X/5,  is  dependent  upon  the  Reynolds 
number,  with  the  value  of  the  ratio  increasing  dramatically  for  low.. 
Reynolds  number  flows.  The  IKET-computed  values  of  the  outer-region, 
length  scale  ratios  for  the  low  Mach  number  hypersonic  flows  suggest 
that  the  ratio  may  also  be  a function  of  the  Mach  number.  In  their 
discussion  of  the  outer-region  length  scale  ratio,  Bushnell,  Cary,  and 
Harris  (Ref.  9)  ascribe  to  this  ratio  a value  Independent  of  the  Mach 
number,  citing  the  results  of  Maise  and  McDonald  (Ref.  71).  The  data 
examined  by  Maise  and  McDonald  covered  a Mach  number  range  from  zero  to 
five  and  confirmed  Morkovin's  hypothesis  concerning  the  structure  of. 
compressible  turbulent  boundary  layers  for  this  Mach  number  range.  '' 
These  data  say  nothing  about  the  behavior  of  the  outer-region  length  - 
scale  ratio  at  Mach  numbers  above  five.  However,  the  data  of  Horstman 
and  Owen  (Ref.  72)  taken  on  the  cylindrical  portion  of  a cone-ogive- 
cylinder  at  a free-streain  Mach  number  of '7.2  exhibit  values  of  the 
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outer-region  length  scale  ratio  in  the  range  from  0,06  to  0.07.  Thus  ,• 
the  behavior  of  this  ratio  at  high  Mach  numbers  is  an  open  question. 

From  the  results  of  this  analysis  it  is  not  possible  to  ascertain 
explicitly  whether  this  ratio  is  decreasing  in  the  low-hypersonic  regime 
or  whether  the  decrease  is  fortultuously  being  caused  by  the  assumption 
that  all  structural  constants  are  invariant  with  Mach  number.  The 
results  of  the  IKET  analysis  as  well  as  the  experimental  data  of  Horst- 
man  and  Owen  (Ref.  72)  do  indicate  that  such  a dependence  upon  Mach 
number  in  the  low-hypersonic  regions  may  be  present  and  significant. 

More  experimental  work  on  turbulent  boundary-layer  structure  is  required 
to  settle  the  question  of  the  behavior  of  the  outer- region  length  scale 
at  elevated  Mach  numbers. 

The  roughness  element  geometry  produced  by  sand  blasting  the  wind- 

tunnel  model  is  of  a very  different  nature  from  the  sand  grain,  roughness 

for  which  the  turbulence  model  is  postulated.  The  results  of  Dir ling  ■ 

(Ref.  36)  as  presented  in  Fig.  3 provide  a correlation  between. any 

roughness  [with  some  restrictions  as  suggested  by  Reda,  et  al..(Ref.  37)) 

and  sand  grain  roughness.  The  sand-blasted  roughness  observed  via"; 

enlarged  photographs  of  sample  cross  sections  is  random..  Because:of 

this  randomness  the  parameters  needed  to  utilize  the  correlation  in  Fig; 

3 were  at  least  averages  and  only  approximate.  Nonetheless,  the  correr 

lation  of  Fig.  3 provided  an  estimate  of  the  ef  f active j sand  graln-jrough- 

ness  size.  The  Dlrllng  correlation  showed  the  effective  sand  grain 

size,  k , to  measured  profilometer  mean  height,  k , to3be  about  0i5  to 
r a . 

1.0.  In  the  IKET  calculations  presented  in  this  section  a value  of  k^  = 

0.0002  ft  was  used,  as  the  measured  k was  3.5  mil  (0.000292  f t),!and , r.nv 

a 

0.0002  was  a reasonable  mean.  The  IKET  calculations  exhibited  little  ! 
sensitivity  to  k^  if  k^  was  of  sufficient  magnitude  to. result  in  fully  . 
rough  flow  initially.  It  was  also  found  that  transition  could  be: 
Induced  by  roughness  effects  for  x/R^  < 0.25  and  that,  whether  induced 
by  roughness  or  tripped,  the  behaviors  of  the  heat-transfer  rate  for 
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x/R^  > 0.25  were  essentially  the  same.  This  calculation  Illustrates 
some  of  the  extent  of  the  physics  embodied  In  the  IKET  extended  mixing- 
length  analysis; 

3.8  COMPRESSIBLE  TURBULENT  BOUNDARY  LAYERS  WITH 
TRANSPIRATION 

f 

Ablation,  which  often  occurs  on  hypersonic  vehicles,  especially 

during  reentry,  Is  usually  simulated  In  ground  test  facilities  by 

transpiration  (l.e.,  blowing).  Thus  an  assessment  of  the  ability  of 

this >' analysis  to  accurately  predict  transpired  flows  Is  In  order.  The 

form  of  the  IKET  equation  as  given  by  Eq.  (16)  can  be  used  to  determine 

• » 

the  changes,  if  any,  that  must  be  made  to  Eq.  (37),  the  IKET  equation 
used  In  this  analysis,  to  account  for  transpiration  effects.  Only  the 

r 

quantities  evaluated  at  the  surface  (y  ° 0)  must  be  examined,  as  the 
integral  terms  explicitly  contain  v In  an  acceptable^  consistent  form  as 
do  the  edge  (y  = (S)  expressions.  The  expression  q^  appears  as  a multi- 
plier in'i every  one  of  the  quantities  to  be  examined  except  p'V'  and 
— 2 ~2 

(v")  .3>iAt  the  surface  q and  v'  are  identically  zero  since  turbulent 
motions  Is  completely  damped  at  the  wall.  Thus  transpiration  adds  no 
terms' tocthe  IKET  equation  as  given  by  Eq.  (37)  - as  long  as  the 
transpiration- Induced iveloclty  perturbation,  v^.  Is  taken  to  be  Iden- 
tically zero.‘;  Thus-’ the  additional  assumption  that  the  transpiration 
f low :>ls!  laminar ~ at5  the  wall  is  herein  invoked,  1.  e. , v'  = 0. 
rii;  ■ •;  f ' ■ 

The  data  of -Squire' (Ref s.  73  and  74),  taken  at  nominal  Mach  numbers 
of  1.8,  2.5,  and'>3'. 6v  were  used  as  the  basis  for  examining  the  IKET 
analys Is ( tfansplfatibn  results.  These  data  were  taken  in  a blowdown 
supersonic 'wind  tunnel,  with  a sintered  bronze  plate  used  as  the  porous 
wall.  Velocity  profiles  were  measured  at  a number  of  stations,  and 
skin-friction  coefficients  were  deduced  by  means  of  the  momentum  Inte- 
gral 'equation  and  the  measured  profiles.  Nominal  test  conditions  for 
the  data  presented  by  Squire  are  ’ given  in  Table  2.  • ..  . 
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Table  2.  Test  Conditions  for  Squire's  Data 


Mach  Number  p , psl  T , "R 

*^0,“’  o,*’ 

Re  /ft 

oo 

(pv)  / (pU) 

W e 

1.8 

26.7 

’ 527.4 

8 X 10^ 

0.0013 

2.5 

44.7 

527.4 

9.7  X 10^ 

0.0013 

3.6 

119.7 

529.2 

15.1  X 10^ 

0. 00065 

The  IKET  analysis  was  applied  by  placing  the  virtual  origin  for  the 
calculations  to  match  the  values  of  the  momentum  and  displacement  thick- 
nesses measured  for  the  untranspired  case  at  the  first  station  where 
profiles  were  reported  by  Squire.  Once  the  location  of  the  virtual 
origin  for  a given  Mach  nximber  was  obtained,  It  was  considered  Invariant 
with  the  blowing  rate: 


B 


(PU)^ 


(77) 


Presented  In  Fig.  31  for  nominal  Mach  numbers  of  1.8  and  2.5  are 
typical  measured  velocity  profiles  (Squire),  velocity  profiles  computed 
using  the  IKET  analysis,  and  velocity  profiles  as  calculated  by  Fletcher 
(Ref.  75)  using  a conventional  mixing-length  approach.".’  Agreement  between 
the  measured  and  IKET-computed  profiles  Is  satisfactory.  The  calcula- 

tlons  of  Fletcher  show  some  discrepancies  when  compared  with  the  experl- 

^ . .o'l!  ■ 

mental  profiles.  Calculations  were  also  made  using  the  computer  code 
developed  for  this  analysis  operating  In  a conventional  mixing- length 
mode  (i.e.,  A/6  forced  to  be  0.09).  Typical  points  from  this  calcula- 
tion are  shown  by  an  asterisk  (*)  In  Fig.  31.  The  points  fair  directly 
upon  the  results  presented  by  Fletcher,  thus  Indicating  the  importance 
of  the  outer-region  length  scales  for  transpired  flows.  Figure  32, 

y 2 

which  presents  /C^/2a^  Y/6  and  B as  functions  of  distance  down  the  plate 
Illustrates  the  effect  of  using  an  extended  mixing-length  fomulatlon. 

In  the  region  where  transpiration  Is  occurlng,  the  outer-region  length 


Lil  ■i 
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scale  decreases  monotonlcally  from  the  expected  value  of  0.09  to  about 
0.08.  The  velocity  profiles  (as  shown  In  Fig.  31)  generated  with  the 
IKET  analysis  confirm  the  effects  of  this  decrease. 


Figure  31.  Transpired  veiocity  profiles  of  Squire. 


The  skln-frlctlon  coefficients  as  deduced  by  Squire,  calculated  by 
the  present  method,  and  computed  by  Fletcher  are  presented  In  Table  3. 


Table  3.  Skin  Friction  Coefficients  for  Squire's  Data 


H 

00 

C^-IKET 

C^-Pletcher* 

Cj-Squlre'* 

B 

1,8 

0.001408 

0.00126 

0.00140 

0.00130 

2.5 

0.001057 

0. 00092 

0.00100 

0.00130 

3.6 

0.00096 

— 

0. 00099 

0.  00065 

®Ref.  75 

'^Refs.  73  and  74 
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Table  3 indicates; that. the.  IKET  analysis  bverpredicts  i by  a. slight- amount 
the  values.' Indicated  ; by;  Squire  and  that  Fletcher ; consistently  underprer. 
diets  the  values  Indicated  by  Squire.  Actually, , both,  the . IKET, predic-.,. 
tlons  and  the  Fletcher  predictions  of  skin  friction. are. within  the 
accuracy  quoted  by.  Squire.  Nevertheless,  the  IKET-predlcted  skin  fric- 
tion coefficients  deviate  less  from  the  experimental  values  than  do  the  < 
values  given  by  Fletcher.  I ^ 

a 0015 

a 0010 

B 

aooos 

aoooo 
a 10 

ao8 

[f]lKET 

ao4 
ao2 
0, 


4.0  CONCLUSIONS 

^1.'  i ’t  ’ . . 

The  IKET- based  extended-mixing- length  hypothesis  as  developed  and 
assessed  in  this  report  has  been  demonstrated  to  be  an  effective 
technique  for  the  calculation  of  many  different  compressible  turbulent 
boundary-layer  flows.  The  ability  to  calculate  the  outer-region  length 


— 

— 

as  1.0  1.5  2,0 


Figure  32.  B and  [X/$]  ik  et  f®*' Squire's  ,, 

= 2.5  case.  ' 

■•A'  A’- 
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scale  for -turbulent  shear  stress,  has  been  shown  to  be 'important  for  • . 
flows  with  significant  pressure  gradients,  roughness,  and  transpiration, 
and  necessary  for  flows  undergoing  transition  and  relamlnarlzation.  The 
baseline  turbulence  model  employed  and  the  modifications  made  to  the  • 
baseline  model  embody  a great  deal  of  the  "physics"  of  compressible 
turbulent  boundary-layer  flows.  For  situations  in  which  the  applica- 
bility of  conventional  mixing-length  analyses  Is  questionable  or  doubtful 
the  extended  mixing-length  approach  delineated  in  this  report  provides  a 
next  logical  degree  of  sophistication.  ' 
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^ APPENDIX  A 

' DETERMINATION  OF  BASELINE  TURBULENCE  MODEL  CONSTANTS 

Classical  turbulent  boundary-layer  theory^  as  well  as  recent 
experimental  data,  Is  used  to  deduce  plausible  values  of  the  empirical 
constants  In  the  turbulence  model.  The  structural  scales  a^  (the  fully 
turbulent  value  of  the  structural  scale  a^,  a2»  and  a^  are  taken  from 
Bradshaw, -Ferrlss,  and  Atwell  (Ref.  18)  and  Mellorand  Herring  (Ref. 

76): 

“ 0.150 
a2  “ 0.566 
33  “ 0.150 

Rose  and  Murphy  (Ref.  26)  concur  with  the  value  of  Bradshaw  et  al.  for 
a^  and  further  suggest  that  the  probable  validity  of  >k>rkovln's  hypo- 
thesis is  the  reason  why  values  of  the  structural  scales  based  upon 
incompressible  flow  data  may  be  used  for  compressible  flows.  It  is 
further  assumed  that  the  a^  and  a^  scales  remain  constant  across  the 
entire  boundary  layer. 

The  empirical  constant  Is  deduced  by  equating  the  classical 
mixing-length  expression  for  In  the  Inner  region  but  outside  of  the 
viscous  sublayer, 

Mj.  = ? y^  3u/ay  (A-4) 

[where  k Is  the  von  Karman  constant  taken  to  have  In  the  baseline  turbu- 
lence model  the  value  of  0.435  (Ref.  10)]  with  the  equivalent  expression 
for  the  Wolfshteln  (Ref.  24)  model  used  In  this  paper: 


(A-1) 

(A-2) 

(A-3) 
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(where  the  structural  scale  a^  has  been  set  equal  to  the  fully  turbulent 
value  since  Eq.  (A-5)  Is  applicable  outside  the  viscous  sublayer)  and 
thus,  using  the  previously  defined  values  of  a^  and  k,  one  finds  that 

Cy  = = 0.2383  (A-6) 

The  empirical  constant  can  be  evaluated  by  neglecting  all  terms 
In  the  turbulent  kinetic  energy  equation  except  production  and  dissipa- 
tion. Turbulence  production  equaling  turbulent  dissipation  Is  appli- 
cable only  outside  the  viscous  sublayer,  where  it  Is  logical  to  assume  ' 

and  also  are  equal.  This  Is  in  essence  a local  equilibrium  model 
of  turbulence  and  represents  the  well-known  mixing-length  approach. 

From  Eq.  (13)  with  turbulence  production  equal  to  turbulence  dissipation 
[see  section  IV  of  chapter  II  In  the  report  by  Laster  (Ref.  20)  , 


3u/3y  ope 


(A-7) 


which  upon  substitution  of  Eqs.  (35)  and  (36)  reduces  to 


0 < - (Q  % I7) 


(A-8) 


and  which  yields 


Cjj  = (2a^)2/c  = 0.3777  ’ ‘‘^(A-Si) 


♦ c'i.'  rC'L  . 


which  is  In  reasonable  agreement  with  the  results  of  Warsl  and  Mertaugh 
(Ref.  77),  who  suggest  a value  of  0.4525. 


The  behavior  of  the  inner  region  of  the  turbulent  boundary  layer  Is 

dependent  upon  the  damping  functions  [Eqs.  (41),  (42),'’^and  (43)]  of  the 

length  scales  and  the  structural  scale  a^  as  given  by  Eqs.  (38),  (39), 

and  (40).  The  constants  A , A-,  and  A determine  the  behavior  of  the 

u & 

Inner  region  damping. 
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The  numerical  value  of  A 
(Ref.  24),  namely 


was  taken  as  recommended  by  Wolfshteln 


Ay  » 0.0160  (A-10) 

The  value  of  A is  dependent  upon  the  near-wall  behavior  of  1 and  the 

u y 

dissipation  of  turbulent  energy  as  discussed  by  Hassid  and  Poreh  (Ref. 

78). 


The  value  of  the  empirical  constant  A^  Is  deduced  using  the  limiting 
near-wall  form  of  the  turbulent  kinetic  energy  equation  given  by  Jones 
and  Launder  (Refs.  79  and  80), 

3/9y  [vT  (3k/3y)]  - ^ e = 0 (A-11) 

Using  the  value  of  p e given  in  Eq.  (24)  and  writing  for  small 
values  of  y as 

p kl/2  „2 

»»  Ajj  Rt  y = Ap  ^ 2) 

y 


yields,  when  substituted  in  Eq.  (A-11), 


y (3\/3y^)  + (3y/3y)  (ak/9y)  = (Cjj/Ajj)  (yk/y^)  (A-13) 

According  to  Jones  and  Launder,  the  turbulent  kinetic  energy  In  the 
near-wall  region  is  quadratic  In  y;  hence, 


k 


2 


ay  . 


a constant 


(A-14) 


and  Eq.  (A-11)  becomes 

2ay  + 2ay  (3y/3y)  = (ayy^/y^)  (A-15) 
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Eqviatlon  (A- 15)  when  evaluated  at  the  wall  results  In 

Ajj  = Cj,/2  = 0.18885 


(A- 16) 


The  evaluation  of  the  empirical  constant  A Is  begun  by  writing  the 

& 

local  turbulence  Reynolds  number  as 


R. 


P»^T^^7(2aj  ^ y y"*" 


(A-1  7) 


where 


+ p y 


(A-18) 


and  where  Eq.  (25)  has  been  used  to  relate  the  turbulent  kinetic  energy 
to  the  turbulent  shear  stress.  The  local  wall  shear  stress  ls,.denoted 
by  T . Squaring  Eq,  (A-1 7)  and  substituting  Eqs.  (40)  and  (43)  for  a. 

^ I 


gives 


R‘ 


-“'PC-Aa  V] 


(A-1 9) 


which  In  the  near-wall  region  where  Is  small  becomies  '‘^'‘ 

3 (y+)^  (T  /x  ) 

Rf  £_w_ 

“ t 2a,. 


(A-20) 


Using  the  near-wall  approximations 
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in  conjunction  with  the  turbulent  shear  stress  as  defined  by  Eq.  (35) 
results  In 


2A  R a^ 
a t t 


(A-23) 


which  is  valid  In  the  near-wall  region.  Equation  (A-23)  via  use  of  the 
friction  velocity 


u 

w 


(A-24) 


and  y as  given  in  Eq.  (A- 18)  can  be  arranged  as 


't'v  ■ '‘t  (Su+/a/)2  <A-25) 

In  the  near-wall  region,  the  classical  work  of  van  Driest  indicates  that 


Tt/T,,  <2  (y+)^/(A+)2 


(A-26) 


3u'^/3y'''  w 1 - (y‘‘')^/(A‘^)^ 


(A-27) 


where  A Is  a numerical  constant.  Use  of  Eq.  (A-26)  in  Eq.  (A-20)  gives 


, a3 

St  n 2 


(A-28) 


2a,.  (A+)‘ 


whereas  using  Eqs.  (A-26)  and  (A-27)  In  Eq.  (A-25)  results  In 


2a^  A 
t a 


2/  +^2 
K (y^) 


A^  (A+)^  [1  - <2  (y+)V(A+)2]2 


(A-29) 
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The  expression  as  given  In  the  above  equation  allows  Eq.  (A-20)  to  be 

written  In  terms  of  A as 

a 


(a  c a"*^ 

f 

1 - K^(y'^)^(A^)^ 

6 7+6 

" > (A-30) 

2 . +.21 
< (y  ) 

3 

2a^(A'^)^ 

which  when  evaluated  at  the  wall  reduces  to 

A^  = (A+/2a^(c)  (Ay  . o.0469  (A-31) 

where  A Is  taken  to  have  the  value  26.0* 

The  last  constant  to  be  determined  is  The  classical  mixing- 

length  expression  for  the  turbulent  viscosity  In  the  outer  region  is 

6^  Ou/3y)  (A-32) 

where  X Is  usually  taken  as  0.09.  Assuming  local  equilibrium  so  that 
and  are  equal  permits  the  turbulent  viscosity  of  Wolfshteln  to  be 
written  as 


- (Cy  /2aj.)  p X^  6^  Oa/3y)  (A-33) 

Comparison  of  Eqs.  (A-32)  and  (A-33)  results  In 

X - /za  fC^  X - 0.2069  (A-34) 

D t p 

based  on  a value  for  X of  0.09. 
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Damping  constant  for  structural  scale 
Damping  constant  for  dissipation  length  scale 
Roughness  element  windward  surface  area 
Roughness  element  cross-sectional  area 
van  Driest  damping  constant,  26.0 
Effective  damping  constant  for  mixing  length 
Damping  constant  for  turbulence  shear  stress 
Constant  in  Eq.  (54) 

Structural  scale  given  by  Eq.  (25) 

Structural  scale  given  by  Eq.  (26) 

Structural  scale  given  by  Eq.  (27) 

Fully  turbulent  value  of  structural  scale  a^ 


; V 1 ' ■ ( ; 

Blowing  coefficient  given  by  Eq.  (77) 
Constant  in  Eq.  (54) 

Dissipation  constant  used  In  Eq.  (24) 

■ r ' - • 

Drag  coefficient  of  roughness  element 
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C 


C 

P 


D 


pk 

IKET 


i 


Skin  friction 

Turbulent  viscosity  constant  used  in  Eq.  (23) 

Constant  In  Eq.  (54) 

Constant  pressure  specific  heat 
Turbulent  kinetic  energy  dissipation 
Wall  damping  for  structural  scale  a^ 

Wall  damping  for  dissipation  length  scale 
Diameter  of  roughness  element 
Wall  damping  for  turbulent  shear  stress 
Incident  energy  absorbed  by  boundary  layer 

3 

Dimensionless  energy  absorbed,  E/p  U 

Turbulent  kinetic  energy  generation  attributable  to  roughness 
Roughness  model  constant  used  in  Eq.  (63) 

Mean  static  enthalpy  • jq 

Peak  heating  ratio  in  Eq.  (76) 

Integral  kinetic  energy  of  turbulence 
1 for  axlsymmetrlc,  0 for  two-dimensional 
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K 

k 


k 


a 


k 


r 


k 


t 


i 


M 


rassetij- 


Pk 

Pr 


Pr. 


Acceleration  parameter  given  by  Eq.  (74)  ,, 

1/2[(uO^  + (y')^  + (w')^] 

Roughness  element  height,  actual 
Roughness  element  height,  sand  grain  equivalent 
Turbulent  thermal  conductivity 
Tensor  Index  In  Eq.  (14) 

Length  scale  for  dissipation 
Center-to-center  roughness  element  spacing 
Length  scale  for  turbulent  shear  stress 
Mach  number 
Edge  Mach  number 
riguo  .Tensord index j In  Eq,  (14) 

Exponent  in(Eqi  (73) ^ 0.85 
Peak  heating  ratio  in  Eq.  (76) 

Laminar  Prandtl  number, ^0. 71 
Turbulent  Prandtl- number,  0.  90,  n 
Fluctuating  pressure  j-.-ow.-:  roi  ? 
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p Mean  pressure 

p^  Dimensionless  pressure  gradient 

p^(x)  Weighed  pressure  gradient  given  by  Eq,  (55) 

Constant  In  Eq.  (66),  0.1 
q Heat-transfer  rate 

q^  Heat-transfer  rate  at  stagnation  point 


2 

q 

R 


Re 


Re 


'e,x^ 

Re  „ 
e,6 

Re  /In. 


St 


(uO^  + (v")^  + (wO^ 

Gas  constant 

Nose  radius 

Reynolds  number 

Edge  Reynolds  number  at  transition 

Edge  Reynolds  number  based  on  momentum  thickness 

Free-stream  Reynolds  number  per  Inch 

-if  - 

Turbulent  Reynolds  number,  pk  y/y 
Roughness  Reynolds  number,  u^k^/v^ 

Radius  for  axlsymmetrlc  flow 

Stanton  number  based  on  free-stream  conditions  and  (T 


- T ) 

oo  w 
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T 

f 

U 

e 

u 

u 

oo 

u' 

+ 

U 

+ 

u 

w 

V 

V 

v" 

+ 

V 

w 

w" 

X 


*2 
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Temperature 

Mean  temperature  ■ ; 

Outer-edge  velocity  , 

Mean  streamwise  velocity  in  boundary  layer 

Free-stream  velocity 

Fluctuating  streamwise  velocity 

Friction  velocity,  u/u^ 

Wall  friction  velocity, 

Velocity  expression  given  by  Eq.  (5) 

Mean  surface  normal  velocity  in  boundary  layer 
Fluctuating  normal  velocity 
Dimensionless  normal  velocity  at  wall 
Fluctuating  lateral  velocity 

Streamwise  coordinate 

Transition  location 

Initial  point  of  laminar  ..flow  breakdown,' 
p^Initial  locatiqnfpf  turbulent- burst 


97 


AEDC-TR-77-96 


Y 

y 

^cr 

+ 

y 

a 

B 

B 

Y 

6 

6* 

* 

6y 

e 

e 

K 

A 


Outer-region  length  scale  for  turbulent  shear  stress 

Normal  coordinate 

Location  of  maximum  turbulent  shear  stress 
Dimensionless  normal  coordinate  given  by  Eq.  (A- 18) 
Constant  used  in  Eq.  (A- 14) 

A « 

Dimensionless  pressure  gradient,  6 (dp/dx)/T 


Weighed  dimensionless  pressure  gradient  given  by  Eq. 
Ratio  of  specific  heats 


Boundary-layer  thickness 


Boundary- layer  displacement  thickness 
Kinematic  displacement  thickness 


Infinitesimal  length  Ih  coordinate  y 
D/p  given  In  Eq.  (23) 

Variable  of  Integration  In  Eq.  (55) 
Boundary-layer  momentum  thickness 


(59) 


^ > 

von  Karman  constant 


Roughness  element  geometry  parameter. 
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- Conventional  outer-region  parameter  in  Eq. 
Outer-region  empirical  constant  in  Eq.  (46) 
Lag  length  parameter  used  in  Eq.  (55) 

Mean  molecular  viscosity 
Turbulent  viscosity  given  by  Eq.  (22) 

Mean  molecular  kinematic  viscosity 
Mean  density 
Fluctuating  density 
Lag  parameter  defined  by  Eq.  (57) 

Shear  stress  tensor 

Shear  stress 

Expression  given  by  Eq.  (56) 
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rs  V' 

e Outer  edge  of  boundary  layer 

■-  r,'  ■ ' 

o,<»  Free-stream  stagnation 

r Roughness 

t Turbulent 
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W Wall 

eo  Free-stream 

SUPERSCRIPTS 

Fluctuation  quantity 

— Quantity  averaged  with  respect  to  time 
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